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Abstract The heating mechanism at high densities during M-dwarf flares is 
poorly understood. Spectra of M-dwarf flares in the optical and near-ultraviolet 
wavelength regimes have revealed three continuum components during the im¬ 
pulsive phase: 1) an energetically dominant blackbody component with a color 
temperature of T ss 10^ K in the blue-optical, 2) a smaller amount of Balmer 
continuum emission in the near-ultraviolet at A < 3646 A and 3) an apparent 
pseudo-continuum of blended high-order Balmer lines between A = 3646 A and 
A « 3900 A. These properties are not reproduced by models that employ a typical 
solar-type flare heating level of < 10^^ erg cm“^ s“^ in non-thermal electrons, 
and therefore our understanding of these spectra is limited to a phenomenological 
three-component interpretation. We present a new ID radiative-hydrodynamic 
model of an M-dwarf flare from precipitating non-thermal electrons with a 
large energy flux of 10^^ erg cm“^ s“^. The simulation produces bright near¬ 
ultraviolet and optical continuum emission from a dense (n >10^^ cm“^), hot 
(T « 12 000 — 13 500 K) chromospheric condensation. For the first time, the 
observed color temperature and Balmer jump ratio are produced self-consistently 
in a radiative-hydrodynamic flare model. We find that a,T k, 10^ K blackbody¬ 
like continuum component and a small Balmer jump ratio result from optically 
thick Balmer (oo —>■ n = 2) and Paschen recombination (oo n = 3) radiation. 
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and thus the properties of the flux spectrum are caused by blue (A ~ 4300 A) 
light escaping over a larger physical depth range compared to red (A ~ 6700 A) 
and near-ultraviolet (A ~ 3500 A) light. To model the near-ultraviolet pseudo¬ 
continuum previously attributed to overlapping Balmer lines, we include the 
extra Balmer continuum opacity from Landau-Zener transitions that result from 
merged, high order energy levels of hydrogen in a dense, partially ionized atmo¬ 
sphere. This reveals a new diagnostic of ambient charge density in the densest 
regions of the atmosphere that are heated during dMe and solar flares. 


1. Introduction 


It is notoriously difficult to explain the origin of heating at high atmospheric 
densities during solar and stellar flares. The broadband color (white-light) dis¬ 
tribution of the ultraviolet and optical emission during large and small flares on 
chromospherically active M-dw a.rf (dMe) stars exhibits the general shape of a 
T = 8500 — 9500 K blackbody ( Hawlev and Fish^ . 1 19921 : iHawlev et all 2003) 
without an indication of a significant Balmer jump at A = 3646 A. A recent 
homogeneous analysis of flare spectra around the Balmer jump has confirmed 
the presence of a hot blackbody (or blackbody-like) component with a color tem¬ 
perature of T > 10,000 K that contributes most o f the radiated near-ultraviolet 
(NUV) and optical e nergy in the impulsive phase ( Hawlev and PetterserJ . [l9^ : 
Kowalski et al . . 20131 . hereafter, K13). Intriguingly, the presence of hot blackbody 


emission with T ss 9000 K has recently also been inferred from Sun-as-a-star 
obser vations during large (X-class) and small (C-class) solar flares ( Kretzschmarl . 
l201l[ l. The hot (hereafter, “hot” refers to temperatures between 8500 and 15,000 
K) blackbody continuum component is often interpreted as evidence of impulsive 
heating in a region with a high density of rie ~ 10^® cm“^ or more, which 
has bem^^nferre^^onyDhenomenologica^jnodelinE_o£_dMe_^ares_(Cram and 
Woods, 19821 : Houdebina 19921 : Christian et ali 200311 . A recent observation that 
directly supports this idea is the Vega-like flare spectrum det ected during the 
2009 Janua ry 16, megaflare event on the dM4.5e star YZ CMi ( Kowalski et all . 
l201fll . 2013h . which indicates that photospheric-level densities can be produced 
and maintained at sufficiently high temperatures during dMe flares. However, 
heating high densities requires a large amount of volumetric energy deposition 
such that the heating per unit mass is non-negligible. Direct heating of the 
photosphere by non-thermal (NT) deka-keV electrons is problematic because 
the electrons experience rapid energy loss as soon as they impact the top of the 
“thick-target” chromosphere. 

A possible mechanism for increasing the density in t he upper chromospher e 


was proposed in the early gas dynamic simulations of iLivshits et al\ (1198111 . 


These authors showed that a shock from precipitating non-thermal electrons 
(with an energy flux of 10^^ erg cm“^ 


s ^) produces a rapid downflow and 


compression of the material to a very high (> 10^® cm“^) density. The devel¬ 
opment of downward moving “chromospheric condensations” (hereafter, CC) is 
accompanied by an upward moving “evap oration” of chromospheric material 


as it is heated beyond a million degrees ( Fisher. Canfield, and McClvmonh . 
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1985r ). Chromospheric evaporation and CCs are thought to be the fundamental 
processes in solar and stehar flares driving the frequently observe d Neupert effect 
( Dennis and Zarrol . lOQUi Hawlev et ai . 19951 : Gnedel et all\l99^ . Several studies 
have explored the role of explosive CCs as the source of white -light emission in 


solar and dMe flares (Livshits et oZI. 198ll: Gan et al. 199^: Katsova. Boiko 


and Livshits, 1997il . but the formation threshold of CCs is strongly dependent 
on the initial atmospheric structure, the energy distribution (po wer law index , 
low-energy cutoff, and flux) of the precipita ting NT electrons (iFishei . 1989ll . 
and the^ime-degenden^jonizatioi^^nieliuii^ Abbet^^iU^Jlawle^ ^ 1999 ^Allred 


et ai. l2005l . l2006l l. thus necessitating an accurate and realistic description of 
the NT electron distribution, the pre-flare atmosphere structure, and the level 
populations of important cooling agents. 

Recent radiative-hydrodynamic (RHP) simulations of dMe (lAllred et al 


200fih flares with the RADYN code ( Carlsson and Stein . 1995 . 1997!) explored the 


atmospheric response to an actual NT electron spectrum that was obtained from 
the peak of an X-class solar flare inferre d from Reuven Rama ty High Energy Solar 
Spectroseopic Imager (RHESSl) data ( Holman et aZ.! . lioO.lh . A constant, solar- 
type flare heating level of the NT electrons was employed with injected energy 
fluxes of 10^° erg cm“^ s“^ (FIO) and 10^^ erg cm“^ s“^ (FH)- These simulations 
produced explosive mass motions indicative of chromospheric evaporation and a 
CC, but with densities far too low to produce strong, blackbody-like white-light 
emission. Moreover, the prominent spectral component in these simulations was a 
hydrogen recombination spectrum with strong Bahner recombination continuum 
emission and a smaller amount of Paschen recombination continuum emission. 
In these models, the heating at high density was achieved through Balmer and 
Paschen continuum backwarming of the photosphere, which increased in temper¬ 
ature by 1200 K in the Fll run and produced additional amounts of optical/red 
(blackbody-like) emission. For the range of non-ther mal electron energy fluxes 
thus far explored in RHP models (F9-F11; see also iHawlev and Fisheij . [l994 
lAbbett and Hawlevl . 1999l l. a flare spectrum with hot blackbody emission has 
not been produced through the development of a CC, th rough photospheri c 
heating from Balmer and Paschen continuum backwarming ( Allred et a/.l . [200^ . 
or throughDhotos 2 henc_Jieatin£_^om_JGra^^n^^U^^ackwamnn^_(Hawley 
and Fisher, 19921 : Allred et ai . 200fil . see also the extens ive comparison of RH P 
model predictions and recent observations in Section 9 of Kowalski et 

In addition to hot blackbody emission, the recent spectral observations of dMe 
flares also show evidence for a secondary continuum component in the NUV at 
A < 3700 A. A jump in flux at these wavelengths is not accounted for by an 
extrapolation of a single hot blackbody component that can be fit to the redder 
wavelengths at A > 4000 A. The jump in flux relative to the blackbody flux 
is smaller in the peak phase than in the gradual decay phase, when the color 
temperature of the hot blackbody component also decr eases. A comparison to 
the Fll RHP model spectra from lAllred et al\ ( 2006ll led to the conclusion 
that Balmer recombination c ontinuum emission was needed to account for the 
additional flux in the NUV ( Kowalski et al . 20inh . However, high resolution 
{R « 40 000) echelle flare spectra show that the NUV continuum from A < 3700 
A can exhibit the color temperature of a hot blackbody (« 11 000 K) while not 
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exhi biting any signature of a Balmer discontinuity or edge feature (' Fuhrmeister 
et ai, I20n8h that is indicative of Balmer recombination continuum emission. 
Instead, the highest order Balmer lines are broadened and appear to diminish 
as they converge into a “pseudo-continuum” at wavelengths between A = 3646 
A and the bluest identifiable Balmer lin e, which is typically H15 A3712 or H16 
A3704 for both large solar and dMe flares (IDonati-Falchi. Smaldone. and Falcianil . 
1984lHawlev and Pettersen . 1991 : Kowalski et al . 2010l l. 


The broadening and merging of the high order Balmer lines in fl ares is thought 
to be due to the well-known Stark effec t from ambient charges (ISvestkal . 1 19631 : 
IPonati-Falchi. Falciani. and Smaldonel . Il985h . The Balmer edge itself is also 
thought to be broadened, but thermodynamically self-consistent (where the 
internal partition function converges) modeling of this effect was not available 
until the occupational probab ility theory of pressure ionization was developed by 
Hummer and MihalasI (198^. This improved description of pressure ionization 


from Stark broadening at the Balmer edge is currently employ ed in white dwarf 
atmospheric modeling codes ( Tremblay and Bergeronl . 2009f) . but has not yet 
been implemented in models of solar and stellar flares to determine if the theory 
can explain the NUV/blue continuum, and in particular, the lack of a Balmer 
edge. The Stark broadening of the higher order Balmer lines is potentially a 
useful tool for constraining the environment (i.e., charge density) where the 
(blac kbody) continuum emission is produced in stellar flares ( Cram and Woodj . 
Il982ll . 

What atmospheric conditions and atomic processes at high density produce 
a flare radiation field that self-consistently exhibits all the important continuum 
components and line properties in the blue and NUV? First, we must relate 
the observables in blue/NUV spectra to physical parameters of the atmosphere. 
Then, the spectra can be used as powerful tools to constrain the heating mech¬ 
anisms that produce these atmospheric conditions. In this paper, we explore 
the possible role of a very dense CC for producing white-light emission during 
dMe flares, which reveals a meaningful interpretation of the spectra that will be 
used to constrain targeted modeling of individual solar and stellar flares. These 
simulations are timely due to new spectral const raints in the blue/NUV and 


optical wavelength regimes (IKowalski et all 120131 1 and modern computational 


facilities enabling calculation of detailed spectra with sophisticated radiative 
transfer codes. Section [2] describes our method of solution and the starting 
atmospheric conditions. Section [3] describes the atmospheric dynamics and a new 
interpretation of the 10^ K blackbody color temperature. Section |4] describes the 
improved modeling of the Balmer edge. Section 0 explores the viability of the 
model, Section E] summarizes the study and presents several conclusions, and 
Section [7] describes future modeling work. 


2. Method of Solution and Initial Atmosphere 

2.1. NLTE Radiative-Hydrodynamic Modeling with RADYN 


We performed ID radiat i ve-hydrodynamic mod eling using the RADYN code 
( Carlsson and Stein . 1992 . 1994 1995, 1997Ll2002h . RADYN solves the equations 
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of mass, momentum, internal energy, non-equ ilibrium level populations, c harge 
conservation, and radiative transfer (see also Abbett and Hawlev . 199 91. The 
electr on and ion temperatures are equal and an adaptive grid ( Dorfi and Drurvl . 


19871 ) with 250 grid points resolves shocks in the atmosph ere. The NLTE (non- 


LTE) problem is solved in RA DYN using the technique of Scharmer ( 198 if ) and 
Scharmer and CarlssonI 1 1985 1. The NLTE level populations for three atoms are 


calculated: a hydrogen atom (six levels including H ii), a helium atom (five levels 
for He I, three levels for He ii, and He ill), and a singly ionized calcium ion (six 
levels including the Ca ill ground state), giving a total of 22 bound-bound (b-b ) 


and 19 bound-free (b-f) transitions calculated in detail (see I Allred et al . 20051) . 


Continua from elements other than H, He, and Ca ii are t reated as backgro und 
continua in LTE using the Uppsala atmos pheres program dGustafsson . jj*^). A 


radiative cooling function from CHIANTI (iDere et al.l . Il997l : Young et all 2003 ) 
at 42 temperature points between 10^ K and 10® K is included to account for 
optically thin cooling. Radiation from several ionization states {e.g., Mg ii, Fe 
II, and Si ii) that are likely optically thick at 10^ K (Hansteen, 2012, private 
comm unication) may be important in flares dAvrett. Machado, and Kurucj . 
Il98(il ) but is not yet included in the models (see Section [T]). 

Each line profile is calculated with up to 100 frequency points across the 
transition using a Voigt profile and complete redistribution. For Lyman tran¬ 
sitions, the effects of partial redistribution are approximated by truncating the 
line profiles at ten Doppler widths. A change to the line broadeni ng due to the 


Stark effect is implemented in RADYN compared to the mod els oflAllred et al. 
( 20061 ). Here, we use the analytic formulae (Method 1) from SuttonI (|l978l ) to 
add the Stark damping parameter to the total damping parameter in the Voigt 
profiles for hydrogen. The Stark broadening of calcium and helium transitions 
remains the same as in previous versions of the code. 

Backwarming from coronal X-ray and extreme-ultraviolet (XEUV) radiation 
i s trea ted a ccording to the method of Hawlev and Fisher ( 1994 ), Gan and Fang! 
lll99nl) . and I Abbett and Hawlev! (Il999l ). which approximates the radiation source 
using a plane-parallel geometry within 1 Mm and as a point-source beyond 
this distance. The emissivity at 41 temperature points from 10^ — 10® K in 14 
wavelength bins from A = 1 — 2500 A, i s gene rated from the Astrophysical Plasma 
Emission Code ('APEC: ISmith et aZ.1 . 120(H ). The hydrogen, helium, and metal 
cross sections are obtained from CHIANTI. 

We ran flare simulations with the sta rting atmospheric p arameters and bound¬ 
ary conditions similar to those used in Allred et all ( 20061) . in order to facilitate 
the comparison. A preflare, plane-parallel ID atmosphere with log g = 4.75 was 
relaxed to have velocities < 10 m s“^ and converged NLTE populations. The up¬ 
per temperature boundary was initially fixed at 6 x 10® K (rie ~ 2.6x 10^® cm ®), 
and constant non-radiative heating was applied at column mass m greater than 
3.16 g cm“^ (tsooo = 0.06) to emulate a convective energy flux. The atmosphere 
was relaxed using a plane-parallel approximation of the X-ray backwarming, as 
would be appropriate from extended overlying active region loops. The effective 
temperature of the relaxed atmosphere is 3630 K, corresponding to that of an 
early to mid M-dwarf. The height coordinate (z) of the atmosphere lies along 
the length of a semi-circular magnetic loop of half-length L 1/2 = 10® cm (10 
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Mm), with z = 0 at the photosphere (T5000 = 1-0, log m = 1.14 g cm“^). A 
modificat ion to the gravitational a cceleration is incorporated near the top of 
the loop ( Abbett and Hawlev . 1999h : otherwise, the geometry is assumed to be a 
vertical tube of constant cross-sectional area. The imposed ID geometry assumes 
the plasma motion is confined to the direction of the magnetic field, as would be 


appropriate for a low plasma j3 (plasma j3 = 


, but is not necessarily satisfied 


in the photosphere of the preflare atmosphere (the equipartition magnetic field 
strength is T2JcGat^he_2hotos2herej_consistent^itli^bservations_(Johns-Krull 
and Valenti, Ififidli l or in the chromosphere durinsf the flare. Detailed molecular 
transitions that are important for flux redistribution in an M-dwarf photosphere 
have not yet been included in the opacity. 


2.2. Flare Heating Prescription 

We model the atmospheric heating from a double power law, non-thermal (here¬ 
after, NT) electron beam injected at t he top of the atmos phere. In this study, we 
seek to directly connect to the work of Allred et al\ (l200 (Th . w hich used the thick- 
target parameters inferred from RHESSI (iLin et al\ . 2002li hard X-ray spectra 


durin g the pe ak of the 2002 July 2 3 X4.8 GOES class solar flare (|Holman el al. 


200.jl . see also Ireland et all ( 201.'t l'l. A double power law of electron energies is 
parameterized by Si = 3.0 (the power law index from the cutoff energy to the 
break energy), Su = 4.0 (the power law index at energies greater than the break 
energy), Ec = 37 keV (the “cutoff energy” or lowest energy of the NT electron 
spectrum), and Eb = 105 keV (the break energy between Si and Su)- The value 
of Ec derived from observations is an upper limit, and Ec = 37 keV is moderately 
large compared to typically assumed or inferred values on the order of 20 keV. 
Note also that the break in hard X-ray photon spectra (and thus also the inferred 
double power law shape of the NT electron spectru m) can be the result of the 
non-u niform ionization of the ambient atmosphere (ISu. Holman, an d Deniii^. 
2011 ) or from energy losses to the return-current electric field (iHolmanl . 2012li . 
The parameters of the NT electron spectrum are poorly constrained during dMe 
flares because the hard X-ray emissi on is too faint to be observed except during 
the largest flares ( Osten et al . 2010l l. A powerful superflare on the RS CVn HR 
1099 is cons istent with a hard NT electron spectrum with a single power law 
index <5 « 3 (jOsten et all 120071 1. similar to the relatively hard double power law 


spectrum employed in our simulations. 

Free parameters which may not be well- c onstr ained by RHESSI observations 


alone are the energy flux (jKrucker et all [201l|) and the precise time profile 


(duration) of the flare heating in any given magnetic loop. We extend the range 
of heating fluxes from Allred et al. ( 2006l l to also consider those with 10^^ erg 


cm 


-2 


s ^ (F12) and 10^^ erg cm ^ 




(F13), in addition to 10^^ erg cm“^ 
s^" (Fll). The F13 flux is not much larger than the largest values that have 
been derived during solar flares (INeidig. Grosser, and Hrovad Il994ll , and such 
a level has recently been suggested to be relevant to the brightest solar flar e 


kernels ( Krucker et all I2OII ; Milligan et al . , 2014 ; Gritsvk and Somov , 2014[) . 


The duration of the heating i n a single magnetic loo 
in some cases to < 25 — 60 s (IWang. Fang, and Dini 


in solar flares is constrained 


20071 : ISchriiver et aH 120061 1 


SQLA: kowalski_ms_v2.tex; 23 June 2015; 3:56; p. 6 






















































































White-Light Flare Emission from a Chromospheric Condensation 


and in other cases may be several minutes ( Qiu et all . 120101 : Warren , 2006). A 
wavelet analysis of high-time resolution hard X-ray light curves for 647 solar 
flares observed with the Compton Gamma Ray Observatory (CGRO) indicates 
a range_nithe_dm;ation_^^jridivldual_^ursts_Jrom_Od^_s_^o_^0_^Aschwanden 


et al., 1 


. The shortest burst durations from high-tim e resolution cons t raints 

s (Robinson et al. 


observed in the NUV in dMe flare s range from 




to 10 s ( Mathioudakis et al. . 2006h . For our heating profile, we use a constant 
energy flux for a duration of 2.3 s, after which the NT electron energy deposition 
decreases to 0, and the cooling phase is modeled up to 5.1 s for the F13 model 
(and 30 s for the models with lower beam fluxes). Short heating bursts represent 
the heating of a single loop, which can be superposed in sequence to e mulate an 
arcad e development as has been done for soft X-ray data of solar flares (|Warrenl . 
200611 . 

The elastic collisions between NT electrons and the ambient charged and neu¬ 
tral comp onents of the atmo s phere heat the gas through Coul omb interactions. 
Following Hawley and Fisher ( 1994 1. Abbett and Hawlev ( 1999ll . and lAllred et al. 
(|2005l . [^06ll . the NT electron energy deposition rate is included as an extra term 
in the conservation of energy equation. The Coulomb heating rate as a function of 
h eight is modele d here using the analytic, non-relativistic test-particle approach 
of Emsliel ( 1978l l. which was m odified to include the height -dependent ionization 
fraction of the atmosphere in iHawlev and Fisheil (19941. This prescription is 
based on the collisional thick-target model ( Brown Il971 1. but i t also include s 
an approximation for pitch angle scattering of the NT electrons ( Emsliel . Il978ll . 
The inelastic collisions between NT electrons and the ambient neutral hydrogen 
component_ofthe_^;tmosrfiere_cause_directexcitationmri^jonlzatioi^Mott and 


Massey, 1965ll : we follow Kasparova et al. ( 2009l l and add these NT collisional 
rates to the therrnal_^;atesforMi^Rrogen_jjCCording_^o_^he_£rescrl 2 tionjnFang, 
Henoux, and Gan ( 19931 1^1. The results from Dalgarno and Grifhngl ( 1958ll are 
included in these rates to account for secondary ionizations of hydrogen by 
electrons that are ionized by the beam from the ground state of hydrogen. 


3. Properties of the NUV and Optical Continuum Predictions 
3.1. Comparison to Observables 

Continuum predictions from radiative-transfer modeling are constrained with 
two simple observational parameters: 1) the Balmer jump ratio and 2) the color 
temperature of the continuum at wavelengths longer than the Balmer edge. 
These two parameters were referred to as Xflare and Tbb , respectively, in a recent 
homogeneous analysis of 20 flares on dMe stars from K13. A blackbody function 
is fit to the blue window (BW) wavelength ranges (see K13) of the flux spectrum 


comprehensive study of NT ra tes was presented in the calculations of iRicchiazzil lll982ti 
and iRicchiazzi and Canfieldl <ll983fl . who found that the H I n = 1 —> oo and n = 1 —> n = 2 
rates are important, but that the n = 2 —> oo rate was negligible compared to the n = 2 —> oo 
thermal rates for their range of beam fluxes. 
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from A = 4000 — 4800 A (the blue-optical wavelength regime) to obtain Tbb, 
while Xflare is the ratio of the flux averaged over 30 A centered on A = 3615 
A (blueward of the Balmer jump) to the flux averaged over 30 A centered on 
A = 4170 A (redward of the Balmer jump). These wavelength ranges are least 
affected by emission lines while also being near the Balmer jump (see K13). 
We note that the detected continuum emission shortward of the Balmer jump 
in K13 also included wavelengths between A = 3420 A and A = 3600 A. A 
more secure detection and characterizat ion of this continuum emi ssion would 
include shorter wavelengths (such as in iHeinzel and Kleint . 2014 . for a solar 
flare observed from space), but the signal-to-noise decreases at A < 3600 A due 
to the atmospheric cutoff. The distribution of these parameters was found to fall 
between Xflare = 1.5 — 4 and Tbb = 9000 — 14 000 K during the peak phases of 
nearly all of the flares in the sample from K13. The value of Xflare,peak varied 
with flare type, with the most impulsive flares exhibiting a narrow range of 
1.5 — 2.2 and the gradual flares exhibiting larger values. We measure these same 
parameters from the flux spectra for the peak phases of the model continua for 
the range of beam fluxes studied here. The model continua are calculated at a 
fixed number of points across the NUV and optical, and so we interpolate the 
continua to the wavelengths used in the observations. We also compute the ratio 
of the line-integrated, continuum-subtracted Hy emission line flux to the value of 
the continuum flux at A = 4170 A (hereafter, H 7 EW 4170 ; note, this was referred 
to as H 7 /C 417 O in K13). 

The model predictions for Xflare, Tbb, H 7 EW 4170 are given for each NT electron 
flux in Tabled! For the F13 simulation, the table gives the coarse time evolution 
of these parameters for the impulsive and gradual phases. The evolution of the 
flux spectrum for F13 is shown in FigurefllFI. At early times, the F13 flux spectrum 
exhibits a large Balmer jump ratio and a cool color temperature in the blue- 
optical. However, from t = 1.2 s to 2.2 s, the continuum at A > 3646 A becomes 
bluer as the flux at A = 4300 A increases significantly while the flux at A Ri 6690 
A increases only marginally. This results in an increase in the blue-optical color 
temperature, Tbb, of 2000 K. Also, the Balmer jump ratio decreases from 3.3 to 
2 as the flux at A w 3550 A decreases while the flux at A = 4300 increases. The 
continuum at wavelengths shorter than the Balmer edge also becomes bluer over 
this time interval as the flux at A < 3000 A increases while the flux at A ~ 3550 A 
decreases. In Section [3.31 we will explain how these changes in the flux spectrum 
are related to the increasing hydrogen b-f opacities in the atmosphere. 

At t = 2.2 s, the F13 simulation produces a flux spectrum that falls into the 
regime of the observations for the peak phase of dMe flares with a relatively 
small Balmer jump ratio, Xflare = 2.0, and a blue-optical color temperature of 
Tbb = 9300 K. The F13 model NUV flux between A = 3000 — 3646 A also 
exhibits the color temperature of a hot blackbody with T = 11 500 K. Although 
the signal-to-noise of spectra at NUV wavelengths is lower than the data at blue- 
optical wavelengths (due to the atmospheric cutoff), the color temperature in the 
NUV is generally consistent with the value derived at blue-optical wavelengths 


^The flux spectra for the F13 simulation are available upon request. 
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(c/. Figures 8 — 9 of K13). The H 7 EW 4170 also falls within the range of the 
observed values at flare peak (c/. Figures 11 — 12 of K13). Blue-optical emission 
with a color temperature higher than r::! 6000 K has not previously been self- 
consistently reproduced in RHD models of dMe flares (K13). The Fll and F12 
simulations do not produce properties that are consistent with observations of 
the impulsive phase: the Balmer jump ratios are too large (> 7) and the color 
temperatures are too low (< 6000 K). From Table [I] it can be seen that the 
time-evolution of the characteristics of the F13 simulation is consistent with the 
general time-evolution observed in dMe flares: the Balmer jump and H 7 EW 4170 
decrease in the rise phase and increase in the decay phase. The timescale of the 
burst (5 s), however, is much shorter than the timescales of most dMe flares 
(minutes to hours). Adding up sequentially heated magnetic strands as in the 
spatial development of a flare arcade {e.g., Warren . 20061 : Qhi el al . 2010h is 
outside the scope of this work, but we intend to pursue this in a future paper. 


3.2. Dynamics and Energy Balance in the F13 Atmosphere 

The evolution of the F13 flux spectrum in Figure [T] is tied to the dynam¬ 
ical evolution of the atmosphere. We discuss here general properties of the 
hydrodynamics. 

The NT electrons initially deposit their energy over a broad extent, with 
significant volumetric heating from 275 km (log m = —2.1 g cm“^) to the top 
of the chromosphere at 450 km (log m = —4.05 g cm“^), and there is non- 
negligible NT electron energy deposition (300 erg s“^ cm~^) directly in the 
corona. The maximum energy deposition from the NT electrons (8 x 10® erg s“^ 
cm“®) initially occurs at Zo = 390 km (log m = —3.05 g cm“^), which is 99% 
neutral but with a significant electron density Ue = 10^^ cm“®. Most of the NT 
energy is deposited initially into NT hydrogen ionization at Zo- As the electron 
density increases at Zo an increasing fraction of the NT energy goes into raising 
the temperature since heating of electrons occurs more efficiently than heating of 
neutrals (note, the thermalization of electrons, protons, and neutrals is assumed 
to happen instantaneously in these models). 

Hydrogen becomes completely ionized by t = 0.002 s and thermal ionization 
increases over NT ionization as the temperature exceeds 20 000 K. He i then 
becomes the primary cooling agent as the temperature reaches 40 000 K. At 
t = 0.003 s, the temperature approaches 65 000 K and most helium at this 
height is now singly ionized. From t = 0.003 s to t = 0.01 s the He ill fraction 
increases from 2% to 99.6% at Zx = 425 km as the temperature rapidly increases 
to 140 000 K. At t = 0.016 sec, a localized temperature maximum of 180 000 K 
forms at zx, and the thermal pressure generates upward velocities of 4 km s“^ 
and downward velocities of 2 km s“^. These mass motions are the seeds to the 
explosive condensation and evaporation which follow from runaway heating as 
material is evacuated from this region of the atmosphere. 

Meanwhile at t = 0.016 s, between z = 450 and z = 460 km (at the preflare 
transition region heights), there is a strong upward flow of 30 km s“® resulting 
from the increase in temperature to 90 000 K of lower density material (previ¬ 
ously the upper chromospheric material) between the site of maximum beam 


SQLA: kowalski_ms_v2.tex; 23 June 2015; 3:56; p. 9 

















A. F. Kowalski et al. 


heating (425 km) and the transition region (450 km). At t = 0.2 s, the evapo¬ 
ration consists of two fronts: a high density chromospheric front with maximum 
velocity of 125 km s“^ and a low density transition region front with maximum 
velocity of 250 km s“^. The downward front is a high density “chromospheric 
condensation” (CC), with a maximum velocity of 95 km s“^. The maximum 
densities of the chromospheric fronts have increased from nn ~ 2 x 10^^ cm“^ 
(t = 0.2 s) to uh w 5 X 10^"^ cm“^ {t = 0.4 s) and the temperatures have 
decreased from 50 000 K to 37000 K. At the location of zx, the temperature 
has reached nearly 10 MK by t = 0.4 s from the runaway NT electron heating 
as the density has decreased (material is evacuated) from nn = 5 x 10^^ cm“^ 
(t= 0 s) to Uh = 6 X 10^^ cm“^ {t = 0.4 s). At the upper and lower ends of 
the transition from very hot (T « 400 000 K—10 MK) plasma to T « 35 000 K 
plasma, upward and downward propagating shock fronts are located at the very 
steep (nearly discontinuous) jumps in gas pressure, where the He ill fraction 
changes from 100% to < 30%0 

At t = 0.33 s, the timesteps of the simulation become unmanageably small 
(At Ri 10 — 100 ns) due to the precision requirements of RADYN in calculating 
the population levels of He i and He ii in regions of strong temperature gradient, 
especially at the upper evaporation fronts. At t = 0.478 s, we decrease the 
accuracy of the minor level populations, with accurate solutions provided only 
for level population densities that are > 10“® of the total elemental population 
density (initially, the parameter is 10“^^), while the second derivative of the grid 
weights is increased to weight more toward the higher density (CC) front. The 
result is a dramatic increase in the time step size until t = 5.1 s as the low-density 
(upper) chromospheric evaporation front is smoothed over (the number of grid 
points here drastically decreases and the temperature is no longer maintained 
at 35 000 — 50 000 K). The results for the simulation at z > 450 km are thus less 
accurate, but this approximation is required in order to compute ID simulations 
with very steep and dynamic temperature and density gradients. We focus only 
on the radiating material in the CC and below, since this higher density material 
produces most of the optical emission; the model above 450 km should not affect 
the analysis at the lower heights. 

Figure [2] summarizes the dynamical evolution of the lower atmosphere at 
t = 2.2 s, when the maximum white-light emission is present. The top panel 
shows the velocity of the downward moving CC. At this time in the CC, the 
maximum hydrogen density is 7.6 x 10^® cm“^, or 17.4x10“® g cm“®, which is 
over a factor of 12 greater than the hydrogen density at this height {z = 255.4 
km) before the flare. This mass density is comparable to the preflare density 
at z ~ 145 km, which lies between the temperature minimum and the lower 
chromosphere. The evolution of the speed and temperature of the maximum 
density zone in the CC is shown in the bottom panels. By t = 2.2 s, the speed of 
the downflowing material has decreased to R::50 km s“^ as it encounters denser 
material and higher gas pressure lower i n the a t mosph ere, the equation of motion 
for which was derived analytically bv iFisher ( 1989l l. Over the simulation, the 


®The speed of sound in the CC at t = 0.4 s is 30 km s ^ whereas material is traveling at 
speeds of nearly 100 km s“^. 
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temperature of the region of maximum density in the CC decreases steadily; 
this is a result of larger hydrogen emissivity for lower temperatures while in¬ 
creasing the density in the CC gives a larger volumetric cooling rate. Eventually 
the temperature nearly settles between 12 500 and 13 000 K at t = 2.2 s as 
continuous injection of NT electrons deposit energy in the CC. After the NT 
electron energy deposition decreases to zero at t = 2.3 s, the CC continues to 
propagate downward, increasing its density to 1.3 x 10^® cm“^ at t = 4.0 s. The 
temperature drops to 7000 K at this time (when the NT energy input is cut off), 
which is much larger than the preflare temperature of « 5000 K at this height 
(196.5 km). 

At t = 2.2 s, the A Ri 65 — 200 keV NT electrons encounter a thick-target 
“wall” at the CC and produce a maximum energy deposition rate of 5 x 10^ 
erg cm“^ s“^ (compared to a maximum of < 10® erg cm“® s“^ initially) at 
z = 256.9 km; the full width at half maximum (FWHM) of the energy deposition 
rate is only 0.46 km (compared to 75 km initially). The simulation becomes 
computationally problematic shortly after this time and therefore the NT heating 
rate is cut off at 2.3 s. The energy budget around the CC is shown in Figure |31 
which gives the individual contributions to the change in total internal energy 
per unit mass. Positive quantities contribute toward heating, excitation, and 
ionization; negative quantities contribute toward cooling, recombination, and 
de-excitation. The lower end of the CC from Figure [5] at 239.4 km is indicated 
by a vertical long-dashed line in Figure [31 In the CC, the NT electron energy 
deposition is balanced by the thin and thick (detailed) radiative losses, resulting 
in a net decrease in the internal energy (recombination and cooling). At deeper 
layers than the CC, E > 200 keV NT electrons can still penetrate and heat the 
atmosphere. At z Ri 200 km, there is net ionization and heating from both Balmer 
continuum backwarming (red diamonds) and NT electron energy deposition. At 
these depths, radiation from the optically thin loss function and from the detailed 
transitions does not counterbalance these sources of energy. There are also some 
contributions from compressional and viscous work where the CC impacts and 
“accretes” the stationary material below it. 

3.3. Properties of the Contribution Function for Continuum Emission in the 
Flare Atmosphere 

The CC decreases in temperature and increases in density over time (Figure 
HI). By t = 2.2 s the CC has attained temperatures >12 000 K and increased 
in density over 12 times the preflare density, thus producing a large amount of 
radiative cooling from hydrogen and other elements (Figure |3|). In this section, 
we explore the height dependent properties of the emergent intensity from the 
CC at three continuum wavelengths calculated in detail: A = 3550 A in the 
NUV, A = 4300 A in the blue-optical, and A = 6690 A in the red-optical, 
which are indicated in Figure HJ In Section 13.11 we found that the NUV and 
blue-optical F13 model peak spectrum (blue spectrum in Figured]) is described 
by two blackbody spectra with color temperatures of T = 11,500 K between 
A = 3000 — 3646 A and T = 9300 K at A > 4000 A in relative amounts given by 
the value of the Balmer jump ratio (r:: 2). In this section, we use the contribution 
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function to the emergent intensity ( Masainl . Il986 ; Carlsson , 1991) to explain the 
origin and characteristics of the white-light continuum emission in this model 
spectrum. We use the f ollowing form of the contribution function, C/ = dix/dz 
( Leenaarts et al. . 2013l l. for a plane-parallel atmosphere0: 


C, = ^ = ( 1 ) 

dz ^ 

where Xl' is the height-dependent extinction coefficient (in units of cm~^; not 
to be confused with Xflare the Balmer jump ratio), Tj/ is the optical depth 

lO^cm 

(Iz Xvdz), and is the height-dependent source function. More simply, 
this equation is or the plasma emissivity {j^) at a given height and 

frequency multiplied by the attenuation of the radiation (determined by e”"^") 
as it propagates outward in the direction of /r from that height 0 We discuss the 
atomic processes that contribute to the continuum emissivity in Section [3.41 

In Figure m we show the contribution function to the emergent vertical inten¬ 
sity (/r Ri 1) for A = 3550 A, 4300 A, and 6690 A at t = 2.2 s. The temperature 
profile (dashed line) of the flare atmosphere is qualitatively similar to that of 
the quiescent atmosphere, except that the flare chromosphere is hotter and has 
moved to higher column mass. Radiation escapes from two principal regions of 
the flare atmosphere: a stationary component (see Figure [2|) produces a wide 
maximum in C/ at z w 200 km in the mid-/upper- flare chromosphere below 
the CC, while the CC produces a very thin component peaking at z 256.3 km 
with maxima in Cj of 950, 140, and 120 erg cm“^ s“^ A“^ sr“^ cm“^ (extending 
above the top of the figure) for A = 3550, 4300, and 6690 A, respectively. Smaller 
values of z (defined in Section [2) refer to larger physical depths, whereas larger 
values of Az refer to larger physical depth ranges in the atmosphere. We measure 
the FWHM of the narrow component of Cj to compare the physical depth 
ranges (Azfwhm) over which radiation escapes from the CC. The Azfwhm 
is the largest for A = 4300 A: the Azfwhm values are 1 km, 1.75 km, and 1.25 
km for A = 3550, 4300, and 6690 A, respectively. The contribution function 
in the stationary component of the flare atmosphere is also markedly different 
among these three wavelengths. Whereas A = 4300 A has a local maximum at 
200 km, and 20 — 25% of the emergent intensity originates from layers below 
the CC at z < 230 km, the other two wavelengths have very little emergent 
intensity originating from these depths, with A = 3550 A showing the least. The 
arrows on the top of the figure indicate the layers at which r>, = 1 for each 
wavelength. The tx = 1 layer occurs in the CC for A = 3550 and 6690 A, and in 
the stationary flare layers for A = 4300 A. Therefore, the intensity at A = 4300 


^The differential dIx/dz in the definition of the contribution function is the same as the 
differential describing the variation of the specific intensity with height in the radiative transfer 
equation; however, these differentials are not equivalent. 

is often used to describe the emissivity in units of photons cm“^ s“^ sr“^ Hz“^ whereas 
r)i, is used for the emissivity in units of ergs cm“® s“^ sr“^ Hz~^; here, we use ju to describe 
the emissivity in units of ergs cm“® s“^ sr“^ Hz“^. 
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A originates over the largest physical depth range (largest Az) and from the 
lowest atmospheric heights due to the smallest optical depth in the CC. 

At times earlier than 2.2 s, a larger fraction of the emergent intensity at 
A = 3550 A and 6690 A originates from the stationary component of the flare 
atmosphere because the CC has a lower density and optical depth. The evolution 
of the optical depth for A = 4300 A and 3550 A and the evolution of the column 
densities (N) for hydrogen in n = 2 and n = 3 in the CC (where v < —5 km s“^) 
are shown in Figure[S] As the CC becomes denser and attains temperatures that 
excite hydrogen to n = 2 and n = 3, the probability of photoionization by light 
with these wavelengths increases. The values of the column densities at t = 2.2 s 
of Nnin = ^ > 3 X 10^^ cm“^ and Nnin = 3) > 10^^ cm“^ define useful 
requirementqj for large optical depths (T 4300 = 0.8, rasso = 5) in any future flare 
model of CCs with material near temperatures of 12 000 K. 

The evolution of the contribution function results from increasing, wavelength-| 
dependent optical depths at optical/NUV wavelengths and can explain the evo¬ 
lution of the flux spectrum in Figure [T] fSection 13.11) . The very large and con¬ 
tinuously increasing optical depth at A = 3550 A causes a decrease in the flux 
at this wavelength between t = 1.2 s and t = 2.2 s because the radiation can 
only escape from a very thin layer. At t = 2.2 s, the radiation at A = 4300 
A is still relatively optically thin in the CC (r < 0.8), and therefore has a 
higher probability of escape, giving rise to a larger FWHM of the contribution 
function in the CC (Azfwhm) and also a larger fraction that can escape from 
the stationary component of the flare atmosphere. The flux at 4300 A increases 
from t = l.2s to 2.2s in Figured and results from the increasing density but 
a lower optical depth due to the lower photoionization cross-section and smaller 
population density for n = 3 compared to n = 2. The larger optical depth at 
A = 6690 A relative to A = 4300 A is caused by a larger photoionization cross- 
section (but same population density), which leads to a much smaller increase in 
flux at this wavelength in Figure[I]from t = 1.2—2.2 s. The NUV flux significantly 
increases at A < 3000 A and the color temperature at A = 3000 — 3646 A becomes 
larger over this time interval due to similar reasoning for the increase of the flux 
at A = 4300 A and the increase of the color temperature at A > 4000 A. Higher 
signal-to-noise flare spectra at A < 3600 A would help to rigorously test the 
slope changes predicted by these models at NUV wavelengths. 

3.4. Opacity Components in the Optical Depth 

3.4-1. The Optically Thick Hydrogen Recombination Flare Spectrum for the 
FIS model 

In this section, we quantitatively demonstrate how the properties of the flux 
spectrum at < = 2.2 s in Figure [T] primarily arise from spontaneous hydrogen 
recombination continuum emission by reconstructing the contribution function 


®See Section 13.41 there are additional opacities from photoionization of the higher levels of 
hydrogen, hydrogen free-free absorption, and opacity from the H“ ion that can also contribute 
towards the optical depth at these wavelengths. 
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to the emergent intensity from atomic emissi vity processes. T he b-f opacity (x^) 
corrected for stimulated emission is given bv iMihalaa (1978, Equation 7-1): 


Xi. = ^i{[ni - n*e 


( 2 ) 


where <Ti is the photoionization cross-section of an atom with an electron in level 
i for light of frequency n* = ric x (ni/ric)* where ric is the number density 
of protons and is the density ratio in LTE which can be calculated 

using the Saha equation. We calculate the hydrogen b-f and hydrogen free-free 
(f-f, bremsstrahlung) opacity at A = 3550 A, 4300 A, and 6690 A for the flare 
atmosphere at t = 2.2 s. The total hydrogen continuum (b-f-ff-f) opacity values 
for A = 4300 and 3550 A are shown in the top left of Figure [6] The hydrogen b-f 
opacity dominates for both wavelengths, but the hydrogen f-f is non-negligible 
in the CC at optical wavelengths. At A = 4300 A, the Paschen b-f (n = 3 —)• oo) 
continuum opacity is largest, whereas at A = 3550 A the Balmer b-f (n = 2 —)• oo) 
continuum opacity dominates. The b-f and f-f opacities from the H~ ion are not 
shown explicitly in Figure [6] but are included in RADYN and in the contribution 
function (Figure 0]). The H“ opacity is small («3.5%) compared to the hydrogen 
b-f and f-f opacities in the CC but dominates at the lower temperatures (T < 
9000 K) in the flare atmosphere at z < 130 km. Although the electron density 
is high in the CC, Thomson scattering is relatively very small. 

Most of the CC is so dense that the LTE condition S = B (where B is 
the Planck function) applies and the source function is only dependent on local 
quantities. The level populations at t = 2.2 s for n = 2 and n = 3 are equal 
to the LTE populations, and collisional rates are generally large compared to 
radiative rates, for most of the flare atmosphere {z = 100 — 256.6 km) except 
at the surface of the CC. The level population Up is the LTE value through all 
of the CC. The hydrogen b-f and f-f continuum emissivity, j\, is shown in the 
top right panel of Figure [51 For simplicity, only the spontaneous recombination 
emission is shown. The relative rates of induced to spontaneous recombination 
are very small at these wavelengths but can become larger at redder wavelengths. 
The hydrogen emissivity can be approximated as {xv,b-f+Xi^J-f}Bi^^ through 
most of the CC; but to account for departures from LTE at the s urface of the CC , 
we use the NLTE hydrogen b-f emissivity from Equation 7-2 of lMihal^ ( 1978h . 
The ratio of emissivities (Ja=355o/ja=430o) in the CC is about 9.5, as we expect 
for hydrogen recombination (Balmer and Paschen) continuum emission at these 
temperatures. When r is non-negligible, the wavelength dependent attenuation 
becomes important in the contribution function (Equation (UJ). The attenuation 
factor, is shown in the bottom left panel of Figure El and illustrates how 
the A = 3550 A emissivity is preferentially attenuated at z = 180 — 200 km, 
where r 355 o = 20. Here, T 4300 ~ 1 which is large enough to attenuate some 
of the radiation from these heights (compare the shift in peaks at this height 
indicated by vertical dashed lines in the top right and bottom right panels). The 
factor gives rise to the wavelength-dependent changes in the FWHM of Cj 
(A^fwhm) and in the fraction of the emergent intensity that originates from the 
stationary component ISection [3.31 and Figure 01). 
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The contribution function, C/, is shown again in the bottom right panel of Fig- 
ure[n](note logarithmic y-scale as compared to Figure|T]which is on a linear scale). 
Integrating Cj over height gives the emergent intensity at /i = 0.95. The integra¬ 
tion over height reveals that the Balmer jump intensity ratio (/A= 355 o/dA= 430 o, 
similar to Xflare) is ~2.2, a value indicative of the small Balmer jump in emis¬ 
sion in the flux spectrum in Figure [TJ The emergent optical intensity ratio 
I\= 4300 /1 \=6690 of 2.3 is a value consistent with the color temperature of a 
hot (Ri 9000 K) blackbody. We have thus shown that hydrogen recombination 
(and bremsstrahlung) continuum emissivity between 11 000 — 14 000 K integrated 
over heights where optical depths become large (^4300 > 0.5) and are wavelength 
dependent (from the hydrogen cross-section dependence and thermal population 
density variation) can self-consistently generate the color temperature (Tbb) of 
the blue-optical emission and the Balmer jump ratio (xflare) in agreement with 
observations of the peak phases of dMe flares. 

3.4-2. The Optically Thin Hydrogen Reeombination Flare Spectrum for the Fll 
model 


The contribution functions for the Fll model at t = 2.2 s are shown in Figure 
[TKa-b). At this time, the atmosphere has not yet exper ienced explos i ve eva po- 
ration, which is consistent with the Fll simulation from Allred et al\ ( 2006l j. In 
contrast to the F13 simulation, the Ci profiles for A = 3550 A and 4300 A during 
the Fll simulation extend over a large portion of the flare chromosphere. The 
Cl profiles also have second peaks near the preflare photosphere, which is not 
shown here because the important molecular species and opacities in M-dwarf 
photospheres have not yet been included in RADYN (Section [2]). The optical 
depths for 4300 A and 3550 A do not reach values of 1 until heights near the 
pre-flare photosphere; the values of T 3550 in the flare chromosphere range between 
0.001 and 0.1, with r 43 oo about an order of magnitude less. Therefore, the NUV 
and optical radiation from the flare chromosphere is dominated by optically thin 
hydrogen Balmer and Paschen continua emission from recombination due to the 
increased electron density (and temperature). The emergent Fll NUV and opti¬ 
cal spectrum is a combination of mostly unattenuated hydrogen recombination 
emission from z = 250 — 400 km, combined wit h any increased photospheric 
radiation from backwarming ( Allred et al. . 2006lf . The value of Xfla 
(ss 9; Table [T]) and Tbb is low (« 5300 K; Table [T]), 
hydrogen recombination at T ss 10 000 K with negligible wavelength dependent 
attenuation. 

If we use the formul a for spontaneo us therr nal hydrogen recombination emis¬ 
sivity CEauat ion 7-2 of iMihalasl. [19781 see alsolAllen (Il963ll and Wiese. Kelleher, 
and Paquette (Il972ll ) in Equation ([I]), we obtain an approximate formula for the 
ratio of the contribution function at two wavelengths, Ai and A 2 : 


is large 
as would be expected for 


CiiXi) 

Cl{X2) 


Ajnjgbf.i ^ 

Xin'lgbf,2 


( 1 _ 1 ) l.StSXlOVFefAj 




g(r2-Ti)/iJ. 


(3) 


at height z with (Te,ne)^ where g^f are the Gaunt factors and ni and 712 are 
the principal quantum numbers that dominate the opacity in Equation (21) for 
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Al and A 2 , respectivel}!^. When ri ,2 are negligible (as for the Fll in the flare 
chromosphere), Ci{\i = 4300)/C'/(A2 = 6690) for ^ = 1 is « 0.8 at T = 10 000 
K and is not much larger than 1, which causes an optically thin Paschen 

recombination continuum to e xhibit a color temperature of Ri 5000 K ( Kunkej 
1970l : Kerr and Fletcheill2014l . ; Section 9 of K13). When ti ,2 are non-negligible 


(as for the F13 at t = 2.2 s, Section [3.4.111 . the ratio of the attenuation factors 
is > 1 which causes the optical continuum to appear hotter. There¬ 

fore, increasing the optical depth makes the Paschen continuum bluer with 
a higher value of Tbb (apparently hotter), as can be seen by approximating 
the atmosphere as a slab with constant (Te,ne) and integrating over height. 
Note, the correction of the b-f opacity for stimulated emission (second term in 
Equation Q) also affects the color temperature because it decreases the opacity 
at A = 4300 A by a smaller factor than at A = 6690 A. However, the actual 
emergent intensity is an integration over height with varying (Te,ne), which is 
why we model the atmosphere gradients, which are apparent in Figure [21 The 
emergent radiative flux is a Gaussian integration of intensity over five ^ values: 
the intensity at a smaller value of /r gives an even bluer optical continuum and 
smaller Balmer jump ratio due to the factor. 


4. Stark Broadening at the Balmer Edge 


Stark broadening is thought to be an important symmetric broadeni ng mech- 


1963: Wordea 

et al. 

19841: Hawley and Pettersen. 1991: Johns-Krull et all 1997: 

Paulson et al. 

. 2006 

Allred et all 200(31). in addition to Dossible turbulent (Eason 

et al.^ |l992|: 

Dovle et all 198S^ and thermal contributions. The linear (first 


order) Stark effect is the splitting of the degenerate orbital angular momentum 
(Z) states of each principal quantum state (n) of hydrogen, due to a net electric 
microfield from the surrounding distribution of char ges. To date, the most com¬ 
plete t heories of t h e linear Stark ef fect are those of Vidal. Cooper, and Smith 
( 197fll . 197ll. fl9^. ISeat^ (ll99(Th IStehle and .Tacauemotl (1993^. and Kenple 


and Griem ( 1968f l. which take into account the perturbations from the slowly 
moving (quasi-static) protons and from the rapidly moving (dynamic) electrons. 
The Vidal et al. treatment is known as the unified theory because it provides 
an accurate treatment of the electron perturbations from the line core through 
the far wings. Theories of Stark bro adening typica lly only present results for the 
lowest order Balmer transitions, but Lemke (1993) extended the results from the 
Vidal et al. unified theory to n = 22. Therefore, this Vidal et al. unified theory 
(which is considered the best treatment of the non-over lapping far line wings and 


1997 


of the higher order Balmer lines; see th e discussion in iJohns-Krull et al. 
and in IVidal. Cooper, and Smi^ ( 1970(V) is most useful to the interpretation of 
flare spectra in the Balmer jump regio: 


^The value 1.578x10® in the second exponential is 13.598 eV / fcs. 

^Beng'tsonj TjmnidijjMnlj<ep]j^J]li7jjj3royided_£esultsjipJxyHljjisnigJjieJjieorY_(rfKeDDle 

and Griem ^^O^^ fo^n^tonperam^^n^^lec^ui^ensitv^engtsor^^^QG^rivat^communi- 
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In future work (Section [7]), we will implement the iLeml^ 1 1997l l extension of 
the Vidal et al. results (see Section [3. In RADYN, the analytic approximations 
of the Stark damping from ISuttonI (|l978[ l are added to the total damping in the 
Voigt profile (Section (2]) ■ Solar modeling o f infra red hydrogen lines has shown 
that this is adequate (ICarlsson and Ruttenl . Il992 ). 


4.1. Occupational Probability of Level n 


Here, we address the atomic physics necessary to interpret the merging of the 
rapidly converging higher order Balmer lines and the apparent lack of a Balmer 
discontinuity (or edge) at A = 3646 A in flare spectra. 

To model the Stark broadening o f the Balmer edge, we apply the occupa¬ 
tional probability theor;^dfivelOTe^hi Hurnmer_and^ihala^ 1988|}_an^Dappen, 
Anderson, and Mihalas ( 19871 ), which is summarized t horoughly in the recent 
model ing study of cool DA white-dwarf atmospheres of iTremblav and Bergeron 
( 20091 ). This theory provides a thermodynamically self-consistent approach to 
pressure ionization, which allows the internal partition function to converge. 
Unlike procedures that truncate the internal partition function, occupational 
probability theory includes a continuous transition from bound to free states as 
upper levels of hydrogen converge and Stark energy shifts cause them to overlap. 
Electric field ionization experiments with sodium have shown that enhanced 
ionization rates of electrons in level n occur when there is a critical electric 
microfield strength causing the highe st Stark state of n to cross with the lowest 


Stark state of n-1-1 (IPillet ct all . 119841 ). At a level crossing, electrons can undergo 
a Landau-Zener transition from n to n -|- I as the electric field fluctuates with 


time above and then below the critical value (jZeneii . 119321 ). For electric field 
perturbations that cause level crossing of n and n -I- 1, all states > n become 
dissolved or destroyed, and an electron excited to n can undergo a reverse-cascade 
to the continuum through a series of Landau-Zener transitions. 

The occupational probability, ui„, is the probability of an atom being per¬ 
turbed by a net electric microfield that is below the critical field strength that 
would dissolve/destroy level n: 


iCn (charged) = 


Fi(3)dp 


(4) 


where P(/3) is the probability distribution of net electric microfield strength 
(/? ) and 0 r. r\t is the critical net electric microfield that dissolves/destroys level 
n ( Seatonl . Il990l ). Both the high frequency electron collisions and quasi-static 
proton collisions contribut e to level dissolution, but w„ is primarily determined 
from the ambient protons (iHiimmer and Mihalasl . 1 1 9881 : iTremblav and Bergeronl . 
I2OO9I ). The ambient electrons contribute to the cascade plus ionization process 
through their high-frequency collisions with the electrons that are excited to the 


cation) provided higher order line profiles for the study of iJohns-Krull et al\ lll997tl . but these 
results have not been made publicly available to our knowledge. 
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dissolved states and through the mixing of Stark st ates allowing the reverse- 
casca de of Land au-Zener transitions to procee d (see iTremblav and Bergeronl . 
2 OO 9 I I . Following Tremblay and Bergeronl ( 2009^ , we adopt the Q-fit electric mi¬ 
crofield model ( Navfonov et a/.l . ll999 l of the Hooper distribution to represent the 
probability distribution P{P). In our calculations, we use the total occupational 
probability obtained by multiplying Equation (|T]) by the probability of level 
destruction from neutral collisions, which a re m odeled using the same metho d 
described in Tremblay and Bergeronl ( 2009l l and Hummer and MihalasI ( 1988h . 
For the conditions in flare atmospheres, the total occupational probability is 
approximately equal to (charged), given by Equation ([3]). 


4.2. Continuous Opacities at the b-f Edge 


Dappen, Anderson, and Mihalas ( 1987ll developed the method for incorporating 
the occupational probability theory into b-f and b-b opacities for hydrogen to 
account for the effects of Landau-Zener transitions and the subsequent cascade 
to the continuum. 

The Landau-Zener b-b opacity for a Balmer line, Hn, is the nominal b-b 
opacity multiplied by the occupational probability where W 2 =1 for 

the conditions in flare atmospheres. The Landau-Zener b-f opacity is the nominal 
b-f opacity (Equation d^])) extrapolated to wavelengths longer than the standard 
b-f edge (Ao), multiplied by the probability that an atom experiences a critical 
microfield strength or greater, causing the upper level of the transition to be 
dissolved. For any wavelength A > Ao, this probability is obtained by using a 
“pseudo” quantum number n*: 


n 


* 


/ _ 

''nf XZ'^e’ 


(5) 


where E = 13.598 eV = 2.1787102 x 10 erg and Z is the atomic number 
{Z = 1 for hydrogen). The dissolved probability (or dissolved fraetion) is 


D{X) = 1 - 


Wrt 


( 6 ) 


where Wn* is the occu pational probability for level n* and i = 2 for the Balmer 
lines. As discussed bv lDannen. Anderson, and MihalasI ( 1987ll the pseudo levels 
(n*) are simple extrapolation tools between the limits of Wn ~ 1 and << 1 
that are useful for estimating the critical microfield strength (where the highest 
pseudo-stark state of n* and lowest of n* -I- 1 overlap) that can cause level 
dissolution through a cascade of Landau-Zener transitions. Also, note, there are 
several different way s to est i mate the critical field strength, and we use the 
popular method from lSeai^ ( 1991^1 . 

Now accounting for Landau-Zener transitions, wavelengths longer than the 
Balmer edge (Ao = 3646 A) allow there to be a fraction of atoms in n = 2 that 
can be photoionized by A > 3646 A, whereas without Landau-Zener transitions, 
photons with A > 3646 A could photoionize atoms only with electrons in n > 2. 


SOLA: kowalski_ins_v2.tex; 23 June 2015; 3:56; p. 18 
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The total b-f opacity, with and without Landau-Zener transitions, in a repre- 
sentative layer of the CC in the F13 model is shown in Figure IS] Due to the 
high ambient proton density (5.6 x 10^® cm“^), the dissolved fraction is >99% 
(the occupational probability, Wn, is < 1%) for the upper levels that contribute 
opacity at A < 3730 A. The nearly complete level dissolution results in a direct 
extrapolation of the Balmer b-f opacity to these wavelengths. As A increases, the 
fraction of n = 2 hydrogen atoms that can be photoionized decreases, due to the 
larger critical microfield required to cause destruction of the upper level of an 
atomic transition. This leads to a continuous transition from the Balmer contin¬ 
uum opacity into the Paschen continuum opacity, which dominates at A > 4200 
A. The hydrogen f-f opacity (without Landau-Zener transitions) is also shown 
in this figure, illustrating that it is small but non-negligible {xf-f/Xtot ~ 0.2) at 
redder wavelengths. 

Hereafter, we refer to these new opacities as Landau-Zener opacities, and the 
pseud o-continuum resulting from this opacity as the Landau-Zener continuum. 
Note, Dapoen. Anderson, and MihalasI (1987) referred to this as the “dissolved 
level pseudo-continuum.” In Figure|9l we show how the Landau-Zener b-f opacity 
from A = 3646 A to A = 4200 A changes as a function of rip (ss Ue). Self- 
consistent atmospheric parameters were obtained from representative heights 
in the Fll and F13 models, at various times, and covering a range of proton 


densities, Up = 10^'^ — 5 x 10^° cm“'^. As the ambient charge density increases, 
the Landau-Zener b-f opacity maximum shifts from A = 3660 A to A = 3730 A. 
Also shown here are the occupational probabilities, Wnjw^ (open circles, right 
axis, W 2 = 1), that are multiplied by the nominal b-b opacity for each of the 
Balmer b-b transitions. Thus the opacity from the lines is transferred to the 
Landau-Zener continuum. 


4.3. Landau-Zener Opacity Modeling Using the RH Code 


We use the NLTE RH code (jUitenbroekl . l2001r ) to model the Landau-Zener b-f 
and b-b opacities and emissivities between A = 3646 A and A = 5000 A for the 
F13 model at t = 2.2 s. The RH code is well-suited for modeling the Balmer 
jump region dne_to_^he_abiht^^o_calculate_^verla£ 2 iiiS.Jiansitions_(Rybicki 
and Hummer, 1992h with more wavelenEfth points than in the orisfinal dynamic 
(RADYN) simulation. By assuming statistical equilibrium and hydrogen number 
conservation, we can use an arbitrarily large hydrogen atom to model the merging 
of higher order lines near the Balmer b-f edge. Statistical equilibrium is a poor 
assumption during the early phases of the flare simulation, but at t = 2.2 s, the 
resulting error is smalH. As in RADYN, the linear Stark broad ening in RH is 
treated with a Voigt profile and the approximate formulae from ISuttonI (|l978f l 
are used; a prescription for quadratic Stark broadening is also included. 

We obtain the atmospheric parameters (column mass, temperature, electron 
density, and velocity) at t = 2.2 s from the RADYN simulation with a six- 
level hydrogen atom, and we calculate the spectrum with the RH code using 


^We compared the RH and RADYN predictions at t = 2.2 s for the Balmer lines and continua 
that overlap, and we found satisfactory agreement. 
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a 20-level hydrogen atom. In the RH calculation, Ca ii and He are treated as 
background (LTE) elements. At each depth in the atmosphere, the effect of 
Landau-Zener transitions is modeled by multiplying the Balmer b-f cross section 
by the dissolved fraction (Equation ([6])) and the Balmer b-b line profile by the 
occupational probability {wn/w 2 ), as in iTremblav and Bergeron (20^. The 
NLTE opacity and emissivity is then calculated in RH with these modifications. 
We show the result of the RH flux calculation in Figure [TUI compared to several 
observed flare spectra during the mid rise (S#27k2l, late rise (S#28), and peak 
phases (S#31) of a large flare on the dM4.5e star YZ CMi from K13. The 
observed spectra have been normalized to the flux at A = 4700 A to show the 
range of properties at A < 3900 A and the lack of color temperature evolution at 
A > 4000 A (see K13 for the evolution of the flare flux without scaling). Because 
the flare area is unresolved, the continuum brightness is not constrained by 
observations. Thus, the RH model surface-flux spectrum has been scaled by the 
factor (where Rflare is the projected radius of a circular flare area, and 

d is the distance to the star) to match the A = 4700 A flux of each observed flare 
spectrum. An instrumental convolution with a Gaussian of FWHM = 13.5 A has 
also been applied to the model. We are not attempting to model this flare here, 
but we use these spectra to illustrate that the model spectrum characteristics 
are similar to the data at least at some phases of the flare evolution. 

First, the blue-optical color temperature between the Balmer lines at A > 
4000 A is well fit to a blackbody with a temperature of Tbb = 10 400 K, 
which is shown in Figure |TU] as the dashed line. Note that this color temper¬ 
ature is larger by 1000 K compared to that in Table [I] due to the extension 
of the Landau-Zener continuum and Balmer wings between A = 4000 A and 
A = 4060 Although the color temperature is represented by a blackbody 
with T = 10 400 K, the F13 model radiative flux is larger by a factor of 
1.7 compared to the surface flux of a blackbody with this temperature. The 
radiation temperature (rrad,/i=i,A =4300 ~ 11600 K) is also significantly larger 
than the color temperature (Tbb = 9300 — 10 400 K). Hereafter, we refer to the 
hot blackbody component as the hot “blackbody-like” continuum component, 
since we have identified the emission as predominantly due to optically thick 
hydrogen (Paschen) recombination emission in Section [3.4.1l Second, a jump in 
flux occurs above an extrapolation of the blackbody to the bluest wavelengths, 
and this jump in flux is the same in both the F13 model and the observations 
of the late rise phase spectrum. Third, there is no Balmer discontinuity feature 
present in the model spectrum at A = 3646 A: the model calculation produces a 
featureless Landau-Zener continuum connecting the Balmer continuum emission 
at A < 3646 A to the bluest significant Balmer line in emission, HIO A3798. This 
effect originates from the direct extrapolation of Balmer continuum opacity in 
Figure |S] to these wavelengths where the opacity is due to the upper levels of 
hydrogen that have nearly 100% probability of being dissolved. 

The RH flux spectrum with the opacity effects from Landau-Zener transitions 
(and without instrumental convolution) is shown in Figure [Tl] around the Balmer 


^®The S#’s refer to the time-sequential spectrum numbers over each flare analyzed in K13. 
^^We use the same continuum blue window (BW) wavelength regions and weighting as in K13. 
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edge region; no Balmer lines are at all apparent between A = 3646 — 3740 A, 
and Hll is barely identifiable. The RH prediction using a 20-level hydrogen 
atom without the effects from Landau-Zener transitions (gray line, “no L-Z”) 
is shown and definitively illustrates that the overlap/merging of the Balmer 
wings indeed form a pseudo-continuum between the Balmer lines at wavelengths 
shorter than 4100 A, but less than the calculation with the opacity effects from 
Landau-Zener transitions (black line, “L-Z”). The two calculations disagree most 
concerning the appearance and nature of the pseudo-continuum between 3646A 
and 3720 A. In particular, the “no L-Z” calculation exhibits a “blue continuum 
bump” at A ~ 3720 A and a large positive slope at wavelengths between 3646 
A and 3720 A, whereas the “L-Z” calculation predicts a continuum with the 
same slope as the observed Balmer continuum at all A < 3720 A. The “no 
L-Z” calculation also shows that the H10-H19 lines do not merge enough to 
conceal the Balmer discontinuity, which is just as apparent if we apply a large 
instrumental broadening. Using a hydrogen atom with more levels (we also tried 
a 35-level hydrogen atom) does not improve the modeling result: the opacity 
effects from Landau-Zener transitions are nec essary to completely conceal the 
Balme r edge at these densities. As discussed bv IPappen. Anderson, and Mihalas 
( 1987f l. the slope of the pseudo-continuum in the “no L-Z” calculation between 
3646A and the blue continuum bump is the result of the rapidly decreasing oscil¬ 
lator strength density of the higher-order transitions. In the “L-Z” calculation, 
Stark broadening decreases the oscillator strength density of the b-b transitions 
while replacing oscillator strength density with the Landau-Zener continuum. 

In Figure [T^ we show how the contribution function at A = 3870 A (between 
H8 A3889 and H9 A3835) and A = 4300 A compare with and without opacity 
effects from Landau-Zener transitions. As expected, the contribution at A = 4300 
A does not change. For A = 3870 A without Landau-Zener transitions, the 
opacity is dominated by Paschen b-f opacity. This wavelength is more optically 
thin than A = 4300 A. Hence there is more contribution from the stationary 
component of the flare atmosphere at z = 150 — 230 km. In the case with 
Landau-Zener transitions, the Landau-Zener opacity increases the optical depth 
beyond that of the optical depth at 4300 A (the contribution at z = 150 — 230 
km is less than in the calculation without Landau-Zener transitions) and more 
A = 3870 A photons escape from higher in the atmosphere in the CC (at z = 251 
km) due to recombination of electrons from the n* = 8.3 pseudo-level (Equation 
®) to the n = 2 level of hydrogen. 

We can now associate the model flux spectrum in Figure fTOl to several atomic 
emission processes. At A < 3646 A the flux is almost entirely Balmer continuum 
emission from free electrons recombining to n = 2. This part of the spectrum also 
exhibits a color temperature Tbb > 10 000 K, with a peak in the NUV between 
A = 2500 — 2600 A. The Landau-Zener continuum flux between A = 3646 — 3760 
A is bona-fide Balmer continuum recombination radiation resulting from free 
electro ns recombining from the n* pseudo (dis solved) upper levels of hydrogen to 
n = 2 ( Dannen. Anderson, and Mihalas . flQ 8 tII : i.e., undergoing reverse Landau- 
Zener transitions from the continuum to n = n*, at which point they transition 
to n = 2. The flux between the rapidly converging Balmer lines from HIO to He 
is mostly Landau-Zener continuum emission from electrons recombining from 
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the dissolved upper levels to u = 2, but at these wavelengths there is a smaller 
fraction of atoms in the atmosphere that experience a critical field strength, and 
thus the Landau-Zener continuum emission is less than at A = 3646 — 3760 A. 
There is also some contribution from the far wings of the Balmer b-b transi¬ 
tions. The far wings of H(5 and He and the redmost extent of the Landau-Zener 
continuum emission produce additional emission above the nominal Paschen 
recombination spectrum as red as A = 4200 A (compare solid and dotted red 
lines). At A > 4200 A, Paschen continuum emission from recombination of 
free electrons to n = 3 dominates, with a smaller (Ri20%) contribution from 
thermal bremsstrahlung (free-free) emission. Paschen recombination radiation 
also contributes toward the radiation at A < 4200 A but to a much smaller 
degree than the calculation without including the Landau-Zener continuum, as 
can be seen by the attenuation of radiation at z = 160 — 230 km in Figure 
m An interesting observational effect of including the Landau-Zener continuum 
emission is that the measured color temperature of the flux spectrum increases 
from 9300 K (6 level, without Landau-Zener transitions, dotted red line, Table[T]) 
to 10400 K (20 level, with Landau-Zener transitions, thick red line), improving 
further the applicability of the F13 model to the measured Tbb distribution at 
peak for a variety of flares (K13). 

We note that the lower order Balmer lines (Ha-H(5) can be used to constrain 
the models. The CC in the F13 model generates broad profiles, large surface 
fluxes, and a reverse decrement (more Hy flux compared to Ha flux) in the lower 
order Balmer lines. The large widths of the lower order Balmer lines in Figure 
1101 result from Stark broadening. However, improvements to the implementa¬ 
tion of the Stark effect (see Section ED are required for a detailed comparison. 
Although the F13 model produces approximately 18 times the surface flux in 
Hy compared to the Fll model (and approximately 6 times that in the F12 
model), the ratio of Balmer line to continuum values is notably lower than in 
the smaller NT electron flux simulations (Table [T|). At t = 2.2 s, the relative 
amount of Balmer line to continuum emission produced by the CC is consistent 
with typical observed values in the impulsive phase: HyEW 4 i 70 ~ 20 in the F13 
model compared to HyEW 4 i 70 ~ 15 — 30 in the late rise and peak phase spectra 
in Figure 1101 Additional flaring areas that are either impulsively heated with 
smaller NT electron flux or gradually decaying from a large electron flux may be 
required to more accurately reproduce the time evolution of the emission lines 
in these individual spectra. 

4.4. A New Diagnostic of Ambient Charge Density in Chromospheric Flares 

The opacity effects from Landau-Zener transitions result in Balmer recombina¬ 
tion radiation longward of 3646A and explain the NUV/blue “pseudo-continuum”| 
previously attributed to the merging of Stark-broadened Balmer lines. The oc¬ 
cupational probability theory as employed in white-dwarf atmospheric models 
provides a self-consistent quantum mechanical method of modeling the opacity 
effects from cascading Landau-Zener transitions. This new modeling explains 
the lack of an observed Balmer discontinuity (or edge) in flare spectra, while 
still allowing a Balmer jump in flux above an extrapolation of a blackbody 
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that is fit to the blue-optical wavelengths in the Paschen continuum (A > 4000 
A). Accounting for Landau-Zener transitions gives an emergent flux spectrum 
from the F13 model at t = 2.2 s with continuum emission bridging Balmer 
continuum emission at A < 3646 A to the highest-order Balmer line that is 
detectable, which is HIO or Hll in the case of the F13 model. The microscopic 
picture that is reinforced here is that slowly moving thermal electrons are able to 
recombine to n = 2 and produce a photon with wavelength >3646A (a Balmer 
continuum photon) due to the nearly complete dissolution of the upper levels: 
1 — Wnjwi >0.85 for levels n > 11, in a dense plasma with Up = 5 x 10^® cm“^. 
The success of modeling the “pseudo-continuum” emission from A = 3646 — 3760 
A in Figure fin] as a red wavelength extension of Balmer continuum emission is 
consistent with our conclusion that the white-light emission, in this RHD model, 
is dominated by a combination of hydrogen Balmer b-f, hydrogen Paschen b-f, 
and hydrogen f-f emission. 

In solar flares, a “pseudo-continuum” co mponent betw een the Balmer lines 
from^« 3700,^^^870_^^_was_^TS^jioted_^^_lZirirLlll980|)_an^ Zirin and Nei- 
dig ( 198lll . The current, best treatment of the modeling of the Balmer jump 
region and pseudo-continuum wasa£ 2 lied_^o_several_solai_^m;es_bLDonati- 
Falchi, Falciani, and Smaldone ( 1985l l. who used static, isothermal models with 
rie = 10^^'^ — 10^^ ® cm“^ and Tg = 7000 — 10, 000 K. These models produced 
a “blue continuum bump” at A = 3675 A, which shifts to redder wavelengths 
and increases in brightness for higher electron density, and has been discussed in 
light of a new solar flare observation with broad wavelength coverage spectra in 
the NUV ( Kowalski. Cauzzi. and Fletcher . 201^. The pseudo-continuum bump 
at A Ri 3720 A in Figure [TT] produced in our calculation without Landau-Zener 
opacity effects (Section 14.31) is very lik ely analogous to the NUV/blue bump 
featur e at A ~ 3675 A produced in the Donati-Falchi. Falciani. and Smaldonel 
(119851 1 lower density slab models. As noted by these authors, this feature is due 
to t he merging of Stark broadened high order Ba lmer lines. Another feature in 
the IPonati-Falchi. Falciani. and Smaldonj ( 1985 1 models is a continuum “dip” 
at 3646A dueto_Jhe_smeaning_o£tiie_^alrnfir_edge^Donati-Falchi, Falciani, and 
Smaldone (|l985l l truncated the model hydrogen atom at n = 84 and broadened 
the Balmer edge using the Stark width of the highest Balmer line. For n = 84 , 
the Stark damping parameter is 30 A ( Svestkal . Il963l : ILhns-Krull et at . 1997) 
which results in a flattening and smearing of the edge in their models. Since the 
Landau-Zener continuum modeling from occupational probability theory does 
not predict a continuum bump at A « 3720 A or a dip near the Balmer edge 
(for very high proton densities), this modeling may also provide a more self- 
consistent explanation for the range in Bal mer jump char acteristics observed in 
solar flare spectra (as in the compilation bv iNeidig . 1983[l . while also improving 
the diagnostic capabilities of these spectra. 

Intuitively, occupational probability theory provides insight beyond a model 
that considers pseudo-continuum from merged lines as just a superposition of 
b-b o pacities. In the classical framework ( Donati-Falchi. Falciani. and Smaldonel . 
19851 1. the opacity of a b-b transition is spread over a larger spectral range by 


the Stark broadening, and is then added to the opacity of other Stark-broadened 
b-b opacities at a given wavelength. In the occupational probability framework, 
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the Stark-broadened b-b opacity is replaced by a larger b-f opacity at a given 
wavelength when hydrogen atoms experience a critical microfield perturbation. 
The larger b-f opacity results from the fact that free states are not degenerate 
(in orbital angular momentum), and therefore the oscillator strength density of 
the continuum is not re-distributed over a superposition of Stark states as for 
b-b transitions. This physical difference explains the large discrepancy between 
the two mod el spectra in Figure ITT] from A = 3646 — 3730 A (see also the 
discussion in IPapnen, Anderson, and Mihail. 1987tl. I n sum mary, the opac¬ 
ity sources from IDonati-Falchi. Falciani. and Smaldond ( 19851 1 are continuum 
opacity from a Stark-broadened b-f edge and a superposition of Stark-broadened 
b-b opacities whereas in occupational probability theory, the opacity sources 
are continuum opacity from Stark-dissolved levels and a superposition of the 
remaining undissolved Stark-broadened b-b opacities. 

Modeling the opacity effects from Landau-Zener transitions provides a new 
physical interpretation of several dMe flare spectra in the literature. We’ve al¬ 
ready shown how this modeling with a very large density (si 5 x 10^® cm“^) 
can account for the continuum and line properties in a flare with HIO or Hll 
as the bluest identifiable Balmer line iFigur e [TT)ll. In the continuum pea k (1:15 
UT) spectrum from the dMe flare studied in iGarcia- Alvarez et all 1 2002l l. H9 or 
HIO also appear as the bluest Ba lmer line; however, their gasdynamic modeling 
using the Jevremovic et aLl ( 1991 ) fit to the Balmer decrement indicates densities 
nearly two orders of magnitude lower than the occupational probability modeling 
indicates based on the bluest identifiable Balmer line. The bluest identifiable 
Balmer lines in most dMe flares with moder ate spectral resolution are typ ically 
H15 and H16 both in the im pulsive phase ( Hawley and Pe ttersenl . 199ll l and 
in the gradual decay phase ( Kowalski et a LI. l2010l. l2013h . Even R = 40 000 
echelle spectra of dMe flares sometimes do not show emission in Balmer lines 
bluer than about H15 and instead exhibit a nearly featureles^3 continuum with 
a color temperature of 9 000 — 11000 K at A < 3700 A (Figures 2 and 5 of 


Fuhrmeister et al. . 2008h . It is also notable that o ne peak flare spectrum o f 


a dMe flare showed Balmer lines as blue as H19 (jFuhrmeister et all 120111 1 


which the authors concluded was a result of only a small enhancement of the 
chromospheric electron density in the flare. Since the Landau-Zener continuum 
resulting from occupational probability theory can fully account for similar char¬ 
acteristics of this spectral region in the flare spectra in FigurefTOl we suggest the 
range of characteristics observed around the Balmer jump in stellar flares may 
be explained with the range of Landau-Zener opacities in Figure [51 

Future modeling of particular flares using occupational probability theory 
can be used to infer the magnitude and gradient of the ambient charge (proton) 
density in the densest regions of the flare atmosphere that are significantly heated 
and produce white-light emission. In Section 03] and Figure [51 we showed how 
the Landau-Zener continuum opacity maximum changes and occupational prob¬ 
abilities dissolve the hydrogen Balmer transitions over 2.5 orders of magnitude 
in charge density. We present a simple method using the Balmer decrement to 


Aside from many faint metallic and helium emission lines. 
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estimate the charge density for an optically thick flare atmosphere in LTE (as 
for the F13 atmosphere studied here). We introduce the Landau-Zener Balmer 
decrement, which is the flux in a Balmer line transition (e.^., HIO or Hll which 
are relatively unaffected by strong metallic or helium lines) that experiences a 
large change in the dissolved fraction of the upper level of this transition, divided 
by the flux of a lower order line that is not affected by level dissolution {e.g., 
HyjH. In Table [21 we show how the Landau-Zener Balmer decrement HII/H 7 
changes for different temperatures and over the density range of Up = 10 ^^ 
to rip = 5 X 10^^ cm“^. This Balmer decrement is obtained by multiplying 
Wn=ii/wn =5 (where « 1) by the Balmer decr ement formula for t he op tically 
thick, isothermal, LTE limit in Appendix B of Drake and Ulrichl (1980). The 
decrement sharply drops for rip > 10 ^® cm“^ as n = 11 becomes dissolved (see 
also Figure El). The observed decr ements for HII/H 7 in the two peak spectra 
from iHawlev and Petterseiil ( 1991 1 are 0.26 and 0.22, which suggests a charge 
density of not much larger than 10 ^^ cm“^ but also not smaller than 10 ^"^ cm“^. 
This method could be used as an alternative to the Inglis-Teller relatioiFI. which 
requires data with high spectral resolution. However, since lines and continua are 
formed at different physical depths and depth ranges (as for the F13 atmosphere 
studied here), we plan to improve the accuracy of this spectral diagnostic using 
a grid of NLTE flare atmospheres in future works. 


5. The Viability of the F13 Model 


M-dwarf flares are more energetic than sola r flares (in total). Ty pical white-light 


AYF 

energies of dMe flares are 10^° — 10^^ erg ( Hawlev et al. . 201411 with occasional 


events exceeding 10^"^ erg in the f/-band alone ( Kowalski et all 2010f l. The areal 
extent of flare footpoints cannot be measured direct l^l" , and hard X-r ay emission 
is far too faint except during the largest dMe flares ( Osten et al . 2010f ). Although 
we cannot determine if NT electron fluxes on M-dwarfs are indeed as high as 
10 ^^ erg cm“^ s“^, typical values of « 10 ^^ erg cm“^ s“^ in strong solar flares 
cannot explain the properties of white-light emission during dMe flares, and thus 
we infer that the energy flux is higher. The atmospheric response to an F13 NT 
electron burst provides insight into the radiative transfer processes that could 
explain white-light emission, but is this very high flux a realistic model for dMe 
flares? In this section, we evaluate its viability and conclude that some issues 
still remain. 

First, is an F13 beam even possible during an M-dwarf flare? The energy flux 
available to accelerate particles at the reconnection site is limited by the free 


Assuming that the Balmer line emission does not originate at a different spatial location, 
for example from a larger area with a lower density, than the continuum emission. 

^^See the discussion in (Hummer and Mihail lll988ll about how the theory of occupational 
probability differs from the Inglis-Teller relationship: the Landau-Zener transitions require 
fluctuations of the microfield around the critical field value. 

Emission model fits, e.g., with a blackbody function or RHD model spectra, give an indirect 
estimate of the flare area and energy flux requirements. 
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magnetic energy in the corona, which can be estimated by -^^finflow = 2 x 10^^ 
erg cm“^ s“^, where t'inflow is the rate at which magnetic field flows into the 
reconnection region and may be between 0.01 — 0.1 times the Alfven velocity. 
The energy flux requirement of 2 x 10^^ erg cm“^ s“^ into each magnetic foot- 
point assumes equipartition of energy among protons and electrons and that all 
free magnetic energy goes into accelerating particles. Using the larger estimate 
for inflow velocity, 0.1 x UAifven gives Bcor ~ 1-5 kG as the coronal magnetic 
field (assuming ^inflow = 0.01?;Alfven gives 3.3 kG as the coronal magnetic field) 
necessary to supply this energy flux. Direct (Zeeman broadening) measurements 
of small-scale magnetic fields of M-dwarfs reveal average photospheric magnetic 
flux density of 3 — 4.5 kG ( .Tohns-Krull and Valenti . Il99fil : IPhan-Bao et al 


2 OO 9 I I. but the coronal magnetic field above these starspots has not been di¬ 
rectly measu red. The mo s t wel l-studied active corona is that of the dM3.5e 
star EV Lac l Osten et aLLl200^ . which exhibits a very hot (10 MK), confined 
{L « 10® cm), and dense {ue ~ 10^® cm“®) component in quiescence. For these 
structures, the minimum magnetic field required for confinement is estimated 
to be sal kG. Other values from microwave data of flares and quiescence have 


suggested coronal magnetic fields on M -dwarfs exceeding 1 kG (|Osten et al 


200fil : Smith. Giidcl. and Audard . 2005 1. which have also been inferred from 


magnetic confinement estimates using the Haisch Simplified Analysis of data 
from 


et al.. 120061 1. Thus, the energy requirement for accelerating an F13 NT electron 
flux requires relatively high (kilo-Gauss) coronal magnetic fields but not outside 
the realm of current measurements. 

Although the energy flux requirements may be plausible (given large M- 
dwarf coronal magnetic fields), the electron number flux for the F13 model 
may be problematic due to return-current effects. The electron number flux 
corresponding to our double power law F13 spectrum is 10^° electrons cm~^ 
s“^. A return-current electric field is established from the charge displacement 
of electrons and protons, given that both particles are energized equally in the 
process of magnetic reconnection and retraction. The initial energy lost per NT 
electron from the return-current electric field is 13 keV Mm“^ at the apex of our 


model loop (Equation 29 of iHolmanl . 120121) . The NT electrons with the cutoff 
energy (37 keV) lose all of their energy and are thermalized after penetrating 
only 3 Mm into the atmosphere, after which the flux of ongoi ng NT el e ctrons 
decreases with distance as more NT electrons are thermalized. iHolmanI (120121 1 
derives the evolution of the NT electron flux distribution accounting for losses 
by the return-current electric field. We use these results and find that only 
NT electrons with Finitiai > 85 keV penetrate to the top of the chromosphere 
without thermalizing. Furthermore, the F13 NT electron density is comparable 
to the ambient density (rie = 2.6 x 10^° cm“®) at the apex of our model coronal 
loop; this condition would cause the ambient electrons to be accelerated and 
would result in instabilities and other effects that are not modeled here. 

A more basic issue is that the number flux of electrons available in the corona 
flowing into the reconnection site is only ss 5 x 10^® electrons cm“^ s“^ (using 
rie X O.luAifven, where B = 1.5 kG is obtained from the energy flux), whereas 
the number flux in the F13 beam is 10^® electrons cm“^ s“^. Even if electrons 
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in the coronal loop were globally accelerated, there are only an additional 4 x 
10^® electrons cm“^ available. However, this simple assessment does not 
account for processes that may increase electron density during reconnection 
and retraction {e.g., Longcooe and Gnidonil . 201lh . Furthermore, high densities 
between 10^® — 10^^ cm“® have been inferred for some (non-flaring) coronal 


structures on the Sun (|Feldman et oil 119941 1 and on the flare star EV Lac 


( Osten et al . 2006ll . In summary, the energy and number flux requirements for 
an F13 beam could be met under certain conditions (a dense corona with high 
magnetic field), but the return-current losses would be significant. 

A final issue with the F13 model is the short presence of the hot, blackbody¬ 
like emission, which appears only at the end of the impulsive phase between 
t Rs 1.8 s and t = 2.3 s. We expect a coadded burst spectrum to be more directly 
comparable to spectral observations, which have integration times of 10 s or 
more. The coadded spectrum over the 5.1 s duration of the impulsive heating and 
the decay phases (dot-dashed light blue line in Figure [T]) exhibits a lower color 
temperature (7500 K) and larger Balmer jump ratio (2.7) than the instantaneous 
model spectrum at t = 2.2 s (Table [1}. The coadded F13 spectrum is not as 
representative of the peak phases of the M-dwarf flares as the instantaneous t = 
2.2 s spectrum, which exhibits prominent Tbb ~ 10“^ K blackbody-like emission 
as observed. A longer duration of NT electron energy deposition for 5 — 10 s may 
produce similar coadded and instantaneous spectra, but the simulation becomes 
very computationally intensive past 2.3 s with sustained heating ('Section l3.2ll . 

The coadded 5 s F13 spectrum is more representative of the atmospheric 
conditions in the gradual decay phase and the very early rise phase. The gradual 
decay phase of flares exhibits a Balmer jump ratio of 2.6 — 3.3 and a blue-optical 
color temperature of about 8000 K (K13), which are similar to the observables 
obtained from the coadded F13 spectrum (Table[T]). The coadded F13 spectrum 
is also consistent with a spectral observation from the very early rise phase of a 
flare from K13, which we show in Figure [T5] K13 characterized two additional 
properties of flare spectra that can constrain models of the time-evolution: the 
slope in the NUV at A < 3646 A and the amount of deviation at red wavelengths 
(at A > 5000 A) from a blackbody fit to the blue-optical range at A = 4000—4800 
A. Although observations of both the early rise phase and gradual decay phase 
exhibit a Balmer jump ratio that are larger than during the peak phase, only 
the early rise exhibits a T ^lOOOO K color temperature at A < 3646 A and a 
small deviation from the blackbody at redder wavelengths. These two common 
properties of the average burst model and the early rise phase observation can 
be seen in Figure 1131 However, we expect that flare emission at any particular 
time (even in the very early rise phase) may consist of a superposition of many 
individual bursts. In order to conclude that any individual model can explain the 
time-evolution at any particular time, future modeling work should incorporate 
a treatment of the spatial evolution of flares, such as the development of two 
ribbon arcades and sequentially heated white-light kernels. 
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6. Summary and Conclusions 


We have simulated the dynamic and radiative response of an M-dwarf atmo¬ 
sphere subject to a short burst of a constant, high flux (10^^ erg cm“^ s“^) of 
NT electrons with a moderately high low-energy cutoff (37 keV). This simulation 
produces a heated (T = 12 000—13 500 K), high density (uh > 10^® cm“^) “chro¬ 
mospheric condensation” with Nh(ji = 2) > 3 x 10^^ cm“^, Nniji = 3) > 10^^ 
cm“^, and a total physical depth range of ~ 20 km, resulting in large hydrogen 
b-f and f-f opacities within a short amount of time (2.2 s) after the start of 
non-thermal electron heating. At this time, the “flare photosphere” (T 5000 = 1) 
is located in the CC at a height that is nearly 250 km higher than the pre-flare 
photosphere, and the density of the CC is over an order of magnitude larger than 
the material at this height in the pre-flare. Significant heating to T > 10 000 K 
at high densities (comparable to the density in the very low chromosphere of 
the pre-flare) is achieved and maintained in the CC and also in the stationary, 
slightly less dense regions of the atmosphere R:!50 km below the CC (Figures [21 

SI). 


The atmospheric response in the F13 model provides a new, accurate descrip¬ 
tion of the radiation field in a flare atmosphere with a dense CC, and moreover, 
shows that such an atmosphere can produce bright white-light emission with 
properties that are in good agreement with the observations. The flux spectrum 
from the F13 model at t = 2.2 s exhibits a blue-optical color temperature of 
Tbb ~ 10"^ K between A = 4000—4800 A and a relatively small Balmer jump ratio 
(Xflare ~ 2.0; Table [TJ , both of which are large improvements compared to the 
range of these parameters in lower beam flux models (FIO, Fll, and F12 keeping 
all else the same with the NT electron spectrum). These properties describe 
typical peak spectra of dMe flares, and we showed the striking match of the 
instantaneous F13 model at t = 2.2 s to the observations of one large dMe flare 
with time-resolved spectra in the impulsive phase IFigure fTO]) . The F13 model is 
the first RHD model to reproduce three generally observed characteristics (Tbb , 

Xflare, and H 7 EW 4170 ; Tabled]) self-consistently. _ 

The pioneering gasdvnamic simulations of Livshits et al. (198l[) and Katsova. 
Boiko, and Livshits (|l997[ ) produced a CC with similar characteristics to those 
in the F13 CC, and also At 1 at A = 4500 A and T k, 9000 K, albeit with 
10 — 30 times lower NT electron energy fluxes. It is interesting that at least an 
order of magnitude difference in heating rates gives similar dynamic responses in 
gasdynamic modeling and radiative-hydrodynamic modeling. Our new radiative- 
hydrodynamic models use a more realistic starting dMe atmosphere with a NLTE 
treatment of the entire atmosphere (photosphere, chromosphere, transition re¬ 
gion, and corona), a detailed treatment of the helium level populations (which 
are important for regulating the onset of explosive evaporation; see Section 
and a modern parameterization of the NT electron spectrum from RHESSI 
data. However, a higher cutoff of the NT electron spectrum (as employed in our 
model) is exp ected to require a larger flux in order to produce similar dynamics 
( Fishei] . Il989l l. Finally, the RADYN and RH codes provide a detailed spectrum, 
which we have shown is important for critically testing the models against the 
observations, for revealing new atomic physics in flares around the Balmer jump. 
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for constraining the charge density in the densest atmospheric regions that are 
heated during flares (Sections 14.3114.41) . and for understanding how hydrogen 
recombination appears as a blackbody-like spectrum with T « 10^ K (Section 

321 ). 

By analyzing the contribution function for the continuum intensity, we con¬ 
clude that the values of the slope (Tbb) and Balmer jump ratio (xflare) of the 
emergent optical and NUV continuum intensity in the F13 flare atmosphere are 
due primarily to hydrogen recombination radiation subject to wavelength de¬ 
pendent attenuation from Balmer and Paschen b-f opacities. Laboratory spectra 
of hydrogen at high density (n » 10^^ cm~^) and temperatures of T ss 10"^ 
K ( Wiese. Kelleher. and Paauettel . [1972 1 do not exhibit the blue-optical color 
temperature of a hot blackbody with T Ki 10^ K. The necessary ingredient 
to build up optical depth at this density and temperature, and hence produce 
a r w 10^ K blackbody-like spectrum with a small Balmer jump, is a large 
physical depth range (A 2 ;) of emitting material along the line of sight. Thus, the 
color temperature in flare spectra does not directly relate to the atmospheric 
temperature so much as it is an indication of how the amount of emitting material 
along the line of sight varies as a function of wavelength (indeed, material must 
be heated to T si 10^ K for hydrogen to be significantly populated out of the 
ground state into n = 2 and n = 3 and generate optical depth). In the CC of 
the F13 atmosphere, the physical depth range over which continuum emission 
escapes is Azfwhm ~ 1 — 2 km, which is not present in laboratory experiments. 
At the blue continuum wavelengths with the lowest opacity (e. 5 ., at A = 4300 
A), the physical depth range over which continuum emission escapes is relatively 
large in the CC and extends to lower heights at 2 ; ss 200 km (50 km below the 
CC). The amount of emitting material along a line of sight can also change due 
to variations in the density at the heights where continuum emission escapes; 
in the F13 flare atmosphere, the density is smaller at 2 « 200 km compared to 
the maximum density in the CC (Figure [5]). Thus, the integration of the height- 
dependent properties (as for the contribution function in Figure 0]) is important 
for determining the emergent intensity and flux spectra. 

Anothe r important result is that b y applying the occupational probability 
formulae ( Hummer_and^dihal^ , 19881 : IPappen. Anderson, and Mihalas , 1987 ; 
Tremblay and Bergeronl . 2009l l to NLTE b-f and b-b Balmer opacities, we can 


explain the pseudo-continuum properties in the spectral region redward of the 
Balmer edge (A > 3646 A) and between the Balmer lines from A = 3750 — 3900 A. 
These formulae approximate how wavelengths redward of the Balmer jump ex¬ 
perience Balmer b-f opacity in addition to Paschen b-f opacity (and f-f opacity); 
the atomic processes that allow for this additional Balmer continuum opacity are 
dissolved upper levels of hydrogen from Stark perturbations by ambient protons 
and Landau-Zener transitions of electrons between dissolved levels of hydrogen. 
The amount of Balmer continuum emission redward of the Balmer jump and 
the amount of high order b-b Balmer emission vary inversely, and therefore the 
relative amounts of dissolved and undissolved components in spectra provide a 
new diagnostic of ambient charge density in the white-light emitting, partially 
ionized regions in the lower flare atmosphere ('Section 14.41 Tabled). 
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The heating mechanism at high density in dMe flares is still unresolved 
since the CC required to produce the observed properties of the white-light 
emission results from a NT electron beam that would experience a considerable 
loss of energy from the return-current electric held. The main results of this 
work are therefore the new interpretation of the blue-optical color temperature 
and Balmer jump ratio (as the variations in the physical depth range and op¬ 
tical depth as a function of wavelength) and of the apparent pseudo-continuum 
between the higher order Balmer lines (as Balmer continuum emission due to 
Landau-Zener transitions of electrons between dissolved levels). Combined with 
future modeling efforts using a more plausible heating mechanism (Section [7]), 
these insights will allow the spectral observables from individual flares to directly 
constrain the atmospheric properties in the dense, lower atmosphere. 


7. Future Modeling Improvements 


We plan to develop a grid of short flare bursts using a larger low-energy cutoff 
of the NT electron spectrum, thereby reducing the number flux of NT electrons 
(and thus return-current effects) while allowing a larger fraction of energy to 
penetrate deeper in the atmosphere. A possible alternative is to employ lower 
NT electron fluxes and explore possible re-accel eration mechanisms in the chro¬ 
mosphere in a similar way as done recently by ( Varadv et al. . 20141 1. Including 


additional energy deposition mechanisms , such as NT proton energy deposition 
{e.g., a neutral beam simulation such as in iKarlickv et al\. 20noh or a co mbination 
of energy deposition by other means {e.g., Fletcher and Hudson . 2008ll may also 


be important for producing a similar atmospheric response with a smaller flux 
of NT electrons. In addition, we plan to incorporate the spatial development 
of a flare arcade of sequentially he ated footp o ints ( e.g.. as in the observations 
of the Bastille Day solar flare f rom Qiu_e£^ ([201T ) or the cartoon two-ribbon 
model of the megaflare event in iKowalski et al\ ( 2012h l in order to reproduce the 
observed timescales of dMe flares, which are longer than the 5 s burst modeled 
here. 

In future work with the RH code, we intend to improve the implementation of 
occumtionnl_ 2 ro^abilit 2 Aheor^,and_Starl^_broaRenin^Hubeny, Hummer, and 
Lanz (I 1994 II fully extended the occupational probability formalism to NLTE, and 
we suggest that flare modeling codes developed in the future should similarly 
include occupational prob ability formalism in the non-equilibr ium rate equation 
(following Equation 2.23 in iHubenv. Hummer, and Lanzl . ll994 1 in addition to the 
NLTE b-b and b-f opacity and emissivity. Here, our goal was to demonstrate that 
large improvements in the NUV model spectrum result from implementing the 
simple modifications_tothe,^ 22 £iii^iIi^.£Iiii§ 2 isii.£i 2 §ILilLDappen, Anderson, 
and Mi halas (198ltl andiTrem blay an d BergeronI ( 20^ )3). W e will also follow the 
work of Tremblay and Bergeron] ( 2009ll and Seaton] ( 199nh and renormalize the 
Stark profiles of Vidal. Cooper, and Smith! ( 1973ll using /3crit as an integration 
li mit, which leads to significantly narrower high er order Balmer lines (cf. F i gure 5 
of iTremblav and Bergeron . 2009l l The theory of IVidal. Cooper, and Smith ( 1973l l 
is only valid to first order when the line wings do not overlap, since the Landau- 
Zener continuum accounts for larger perturbations. Therefore, too much emission 
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is probably produced between the Balmer lines between A = 3800 — 4000 A in 
the model spectrum in Figure [TTJ Between A = 3646 — 3800 A, the low values of 
Wn (that multiply into the b-b opacity) lead to only a small amount of blended 
Balmer line wings compared to the Landau-Zener continuum, so the error is 
not important. Finally, the linear S tark broadeni ng theory of the Balmer lines 
using the analytic approximations of SuttonI ( 1978h a s Voigt damping parameters 
should be r eplaced by th e full Stark profiles from IVidal. Cooper, and Smith 
(1973) and Lemkel (1993) the conditions present in flare atmospheres, in 
order to give accurate diagnostics on the electron density for unblended Balmer 
lines. 

In future work with the RADYN code (Allred et al. 2015, in preparation), we 
intend to include the following improvements to the modeled physics: the cooling 
from Mg II and Fe ii ions, an improved optically thin radiative loss function, non- 
thermal ionization of He i and He ii, return-current heating, magnetic mirroring 
of the NT electrons, and the fully relativistic Fokker-Planck solution to energy 
deposition and pitch angle scattering of NT electrons. 
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Table 1. Observables calculated from the model flux spectra. These values 
have been obtained without subtracting the preflare spectrum. 


Model (time step; phase) 

Xfiare 

Tbb (K) 

H7EW4170 

Ell (£ = 2.2s; impulsive peak) 

9.2 

5300 

150 

F12 (t = 2.2s; impulsive peak) 

7.4 

5700 

85 

F13 (£ = 0.2s; impulsive early, rise 1) 

6.0 

6300 

70 

F13 (t = 1.2s; impulsive late, rise 2) 

3.3 

7400 

40 

F13 {t = 2.2s; impulsive peak) 

2.0 

9300 

20 

F13 {t = 4s; gradual decay) 

3.3 

5400 

60 

F13 (average burst) 

2.7 

7500 

30 


Table 2. The Landau-Zener Balmer decrement (HII/H7) in the optically thick, isother¬ 
mal, LTE limit. Here, we give only the charged component of Wn (Equation ((3||, averaged 
for the range of temperatures), ignoring the minor dependence on the destruction of levels 
from collisions with other neutrals. 


Up (cm ^) 

Wn^ll 

T= 10 000 K T = 12 500 K T = 15 000 K r = 20 000 K 

1013 

1.00 

0.30 

0.33 

0.35 

0.37 

loi-i 

0.98 

0.30 

0.32 

0.34 

0.37 

5 X IQi'i 

0.85 

0.26 

0.28 

0.30 

0.32 

1015 

0.69 

0.21 

0.23 

0.24 

0.25 

2.5 X 1015 

0.35 

0.11 

0.12 

0.12 

0.12 

5 X 1015 

0.15 

0.05 

0.05 

0.05 

0.05 
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Wavelength (A) 


Figure 1. The evolution of the NUV and optical flux spectra during the F13 simulation. 
The observable parameters extracted from these spectra are given in Table ^ The coadded 
(average) burst spectrum is discussed in Section [5] The preflare spectrum is not visible on this 
scale. The vertical dotted lines indicate the NUV (A = 3550 A), blue-optical (A = 4300 A), 
and red-optical (A = 6690 A) wavelengths used in the continuum analysis (Section [3]3jl. 
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Figure 2. Top: The velocity field in the lower atmosphere at t = 2.2 s shows a chromospheric 
condensation (CC) with a large hydrogen and electron density exceeding 5 x 10^^ cm“^. Bottom 
left: The time-dependence of the velocity in the highest density zone in the CC. Bottom right: 
The time-dependence of the temperature in the highest density zone in the CC. The velocity 
of the CC changes smoothly through the burst while the temperature refiects the cutoff of the 
NT electron energy deposition at t = 2.3 s. 
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Figure 3. Contributions to the energy balance for t = 2.2 s in the F13 simulation. These 
curves show the net amount for each component. For the detailed radiation, the heating and 
cooling are integrated over the wavelengths of the transitions. The Balmer continuum (red 
diamonds) heating/cooling is shown as a separate curve from the total detailed radiation 
curve (which includes the Balmer continuum). The CC, with velocities lower than 5 km s“^, is 
material between z = 239.4 km (long dashed vertical black line) and 257.4 km (thick vertical 
dash-dotted line). 
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Figure 4. The contribution to the emergent intensity for three continuum wavelengths in the 
NUV, blue, and red at t = 2.2 s compared to the temperature structure (for the same height 
range as in Figure EJ. The column mass scale for t = 2.2 s is shown at the top axis. The 
arrows indicate where r = 1 occurs for each wavelength. The contribution functions in the CC 
(z ^ 239 — 257 km) are extremely thin and extend off the top of the plot (see text). 
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Time (s) 


Figure 5. The evolution of column density for n = 2 and n = 3 of hydrogen (left axis, 
dashed lines) and optical depth (right axis, solid lines) for atmospheric heights in the CC with 
velocities < —5 km s“^. At t = 2.2 s, the CC has an extent of 18 km from 239 km to 257 
km and has temperatures > 12 000 K. Approximately 3.5x10^^ erg cm“^ s“^ of NT electron 
energy flux is being deposited within the CC at this time. 
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Height (km) 


Figure 6. The contribution function (bottom right) for 4300 A and 3550 A is broken up into 
three parts: opacity (top left), emissivity (top right), and the attenuation (bottom left). The 
ratios of hydrogen b-f opacity to total hydrogen (b-f and f-f) opacity are shown in the top left 
panel (right axis). In the stationary flare layers {z 190 km), there is a large emissivity at 
A = 3550 A but also a very large attenuation due to photoionizations of hydrogen from n = 2 
in the higher layers of the CC. The peaks of the emissivity and contribution function in the 
stationary layers for A = 4300 A are indicated by vertical dashed lines (see text). 
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Figure 7. The contribution function for the Fll model at p = 0.95 and t = 2.2 s compared 
to the (a) electron density and temperature and (b) optical depth. Only heights corresponding 
to the flare chromosphere, transition region, and lower corona are shown. 
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Wavelength (A) 


Figure 8. The hydrogen b-f continuum opacity accounting for Landau-Zener (L-Z) transitions, 
at t = 2.2 s in the F13 simulation, at a height of 255.4 km (logio m = —2.57 g cm~^) which is 
the maximum of Cj{\ = 3550)A2: x C/(A = 4300)A2: where As: is the vertical extent of a grid 
cell. The values of Up and Tg are respectively, 5.6 x 10^^ cm“^ and 12,770 K. At this height, 
hydrogen is nearly 75% ionized. Also shown are the hydrogen f-f opacity and the hydrogen b-f 
opacity without accounting for opacity effects from Landau-Zener transitions at this height. 
At this layer, rssso = 1.1 and r43oo = 0.2. 
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Figure 9. The hydrogen b-f opacity accounting for Landau-Zener transitions at representative 
times and heights in the Fll and F13 simulations. Each curve corresponds to approximately 
the ambient charge density at each of the tick marks in the colorbar. The opacity curves 
have been scaled to a common value at A = 3600 A to illustrate the wavelength shift of the 
opacity maximum compared to the Balmer b-f edge (vertical dashed line) without accounting 
for Landau-Zener transitions. The occupational probabilities (right axis) of the upper levels of 
the Balmer transitions are shown as circles (for reference HIO is at A ^ 3800 A). 
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Figure 10. RH calculation with a 20-level hydrogen atom and the b-f and b-b opacities 
including Landau-Zener transitions, during a snapshot of the F13 simulation at t = 2.2 s 
(solid red). A T = 10 400 K blackbody fit to continuum windows at A = 4000 — 4800 A 
is also shown (dashed light blue), along with the RADYN prediction with a 6-level hydrogen 
atom (red dotted line). Instrumental broadening has been applied to the red spectrum. Several 
observations during the rise and peak phases of the IF3 flare of K13 are shown: peak spectrum 
is S#31 from K13 (dark gray), mid rise phase is S#27 from K13 (light gray), and late rise 
phase is S#28 from K13 (black). All spectra are scaled to the flux at A = 4690 — 4710 A. 
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Figure 11. RH spectra of the Balmer jump region without instrumental convolution applied 
and with a magnified fiux scale. Black line: Electric field effects on b-b opacities and Lan- 
dau-Zener transition effects on b-b and b-f opacities from critical electric fields. Hll is the bluest 
identifiable Balmer line. Gray line: Electric field effects on b-b opacities and Landau-Zener 
transitions not included. H13 is the bluest identifiable Balmer line. Including Landau-Zener 
transitions in the b-b and b-f opacities predicts more continuum emission between 3646A and 
3720 A, whereas not accounting for Landau-Zener transitions produces a “blue continuum 
bump” at A Pi 3720 A from the blending of Stark-broadened Balmer lines. Instrumental 
broadening is not shown; it broadens the Balmer lines but does not have a strong effect on the 
amount of continuum between the lines. 
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Figure 12. The contribution function at A = 3870 A and 4300 A accounting for opacity 
effects from Landau-Zener transitions (solid) compared to the original calculation (dashed). 
With L-Z transitions, the additional Balmer opacity at A = 3870 A causes more emission to 
escape from the CC and less from the stationary flare layers than without the L-Z transitions. 
The contribution function at A = 4300 A is not affected. 


SOLA: kowalski_ins_v2.tex; 23 June 2015; 3:56; p. 44 












White-Light Flare Emission from a Chromospheric Condensation 


1.0 


0.8 

X 

D 

E 0.6 

«L) 

^ 0.4 

O 

Z 

0.2 


0.0 

4000 5000 6000 7000 8000 9000 

Wavelength (A) 


Figure 13. The coadded (average) model F13 spectra (from RADYN; instrumental broad¬ 
ening is not applied) over 5 s compared to the early rise spectrum (S#26) during the IF3 
flare from K13. The coadded model spectrum has a larger Balmer jump and lower optical 
color temperature than the instantaneous model spectrum at t = 2.2 s in Figure flOl The flux 
at A < 3646 A increases towards bluer wavelengths in the model and observation. Note that 
the opacity effects from Landau-Zener transitions are not included here because statistical 
equilibrium is a poor approximation during the early time steps in the model (Section [4]3jl. 
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